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SUMMARY 

Previous studies have indicated several benefits, including substan- 
1 tial gust alleviation, for aircraft employing an unconventional wing, free to 
' pivot about a spanwise axis forward of its aerodynamic center and subject only 
to aerodynamic pitching moments imposed by lift and drag forces and a trailing- 
edge control surface. A disadvantage of this basic free-wing concept is the 
relatively low wing lift coefficient available, since there is no mechanism to 
| overcome the powerful negative pitching moments associated with the deflection 
of high-lift trailing-edge flaps. 

The f ree-wing/f ree-trimmer is a NASA-conceived extension of the basic 
free-wing concept intended to provide sufficient trimming power to permit the 
j use of high-lift trailing-edge flaps on free-wing aircraft. The wing is con- 
trolled by the pitching moment produced about its hinge axis by the aerodynamic 
forces on a smaller, external trimmer surface attached to the wing by a boom or 
equivalent structure. Furthermore, the trimmer itself is free to pivot about a 
spanwise hinge axis forward of its aerodynamic center. Pitch control of the 
entire assembly is effected by deflection of a trailing-edge flap on the trim- 
ming surface. 

This report describes an analytical study of the longitudinal behav- 
ior of representative small f ree-wing/f ree-trimmer aircraft. It includes an 
assessment of the response to symmetric vertical turbulence, the nature of the 
characteristic longitudinal stick-fixed modes of motion, and the maximum trimmed 



lift coefficients obtainable. Both forward- and aft-mounted trimming surfaces 
are considered. j 

Although the permissible trimmer moment arm is limited by adverse 

i 

dynamic effects, it is concluded that configurations with trimmer surfaces 
mounted aft of the wing hinge can be designed that provide excellent allevia- 

I 

tion of vertical gust loads yet exceed the maximum lift capability of pure 
free-wing aircraft and meet fundamental criteria for the stick-fixed longitu- 
dinal stability. Forward trimmer configurations were found to have inferior 
gust alleviation characteristics. 

While the allowable center of gravity range for the fuselage assembly 
is comparable to conventional aircraft, the wing/trimmer assembly should be 
mass-balanced about the wing hinge to facilitate smooth landings. Close mass- 
balancing of the free trimmer surface, with respect to its hinge axis, is man- 
datory to preclude short-period dynamic instabilities. 

INTRODUCTION 

Background 

Low wing loading has long been synonymous with poor ride quality in 
turbulence, a fact which has probably been a significant deterrent to wide- 
spread acceptance of light aircraft as a practical means of transportation. 

The problem is compounded by the fact that light aircraft spend a major portion 
of their flight time at the lower altitudes where measurable turbulence is most 
likely to be encountered, even in good weather. 

An increase in wing loading can be expected to produce ride improve- 
ment, but this approach is constrained by the need to maintain low minimum fly- 
ing speeds. Not only is a low speed capability desirable from an operational 
standpoint for aircraft in this class, but Federal airworthiness standards (FAR 
Part 23) require a minimum speed of 61 knots or less for single engine aircraft. 
Even with well designed mechanical high-lift devices, it appears that a wing 
loading of 1900 N/M^ (40 lb/ft^) is an approximate practical upper limit; this 
may be compared to wing loadings of 4300 N/M^ to 5300 N/M^ (90 lb/ft2 to 
110 lb/ft^) typical of the commercial jet transport aircraft which have enjoyed 
wide acceptance by the traveling public. 
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The free-wing concept is a novel and very effective approach to gust 
alleviation which would appear to have other benefits as well, but a disadvan- 
tage of the concept is the relatively low maximum lift coefficient obtainable. 

The f ree-wing/f ree-trimmer is a NASA-conceived extension of the basic 
free-wing concept intended to provide sufficient trimming power to permit the 
use of high-lift trailing-edge flaps on the free wing. 

Free-Wing Concept 

The basic concept of the free wing was disclosed in U.S. Patent 
No. 2,347,230. This patent, now expired, was issued in 1944 to Daniel R. Zuck, 
who built a small prototype aircraft in 1945 as a private venture. This air- 
craft was never successfully flown. 

As conceived by Zuck, a free-wing aircraft differs from a conven- 
tional airplane in that the two panels of the fuselage-mounted wing are free to 
move independently about a spanwise axis and are controlled by means of trail- 
ing-edge control surfaces. Each wing panel is completely free to rotate about 
its spanwise axis, subject to aerodynamic moments but otherwise unrestricted by 
mechanical constraints. To provide static pitching stability, the axis of ro- 
tation is located forward of the chordwise aerodynamic center of the wing panel, 
as shown in figure 1. The wing is brought to an equilibrium angle of attack 
through a balance of moments created by the trailing-edge surface, which is 
controlled by the pilot, and the torques produced by the lift and drag forces. 

The gust alleviation feature of the free wing is caused by the fact 
that a stable lifting surface tends to maintain a prescribed lift coefficient 
by responding to natural pitching moments which accompany changes in flow di- 
rection. While all stable aircraft tend to relieve the lift increment due to a 
vertical gust by pitching into the relative wind, the rapidity of the alleviat- 
ing motion depends upon the pitching moment of inertia. Because of the greatly 
reduced inertia of the wing panel, compared to the aircraft as a whole, the 
free-wing concept produces a significant reduction in turbulence responses. 

Previous work 

The first known analytical study to predict the fundamental dynamic 
behavior of free-wing aircraft permitted independent motion of the left and 
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FIGURE 1. CROSS-SECTIONAL ILLUSTRATION OF THE FREE WING 

right wing panels and is reported in reference 1. The following conclusions 
were drawn from this earlier work: 

(1) Most atmospheric turbulence effects were greatly reduced, par- 
ticuarly in the root-mean-square (rms) load factor (62 percent reduction) and 
rolling disturbances (25 percent reduction). On the other hand, the rms fuse- 
lage pitch rate was increased about 180 percent in comparison with equivalent 
fixed-wing aircraft. 

(2) All stick-fixed modes of motion were stable except for the spi- 
ral mode, where the rate of divergence was found to be excessively high in the 
approach condition. 

(3) The lateral-directional handling qualities were unsatisfactory 
because of the combination of low roll damping and spiral divergence. 

(4) Artificial stability augmentation, in the form of a simple roll 
damper, provided excellent lateral control and turbulence penetration 
characteristics . 

As a result of the very substantial gust alleviation predicted in 
reference 1, a second study was performed to provide a realistic and 
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^comprehensive assessment of the practical aspects of implementation of the 
free-wing concept for light, general aviation aircraft. 

, The second investigation is reported in reference 2. This study, 

I while analytical, was supported by limited wind-tunnel experiments. 

From the results of the second study, the following conclusions were 

j drawn: 

(1) The free-wing concept can be applied to unsophisticated low wing 
loading light aircraft to provide a ride quality, based on normal load factor 
attenuation, equal or superior to that of aircraft employing much higher wing 
loadings. Compared to similar light aircraft in cruise flight, reductions of 
about 54 percent can be realized in the rms load factor increments in contin- 
uous turbulence. 

(2) For free-wing aircraft without differential wing panel freedom, 

I all pertinent handling qualities and certification criteria can be met without 
recourse to stability augmentation, either active or passive. 

(3) Differential pitch freedom between the left and right wing 

; panels should not be permitted for aircraft in this class, although it appears 
1 that passive mechanical devices can be applied to correct the serious lateral 
deficiencies which accompany such freedom. 

■ (4) Leading-edge slats are necessary for takeoff and landing to com- 

pensate for the inherent low maximum trimmed lift coefficients obtained with 
trailing-edge control surfaces. 

(5) The free-wing panels should be balanced about the spanwise hinge 
axis with leading-edge slats retracted. A ballast weight penalty is incurred 
; which might range from 1.5 percent to 7.0 percent of the aircraft gross weight, 
i depending on the detailed design. 

Purpose of this Investigation 

[ The purpose of this work was to perform an analytical evaluation of 

I the dynamic characteristics of an aircraft employing a NASA-conceived concept 
J of a wing and wing-trimmer assembly. In this concept, the wing is free to 
pivot freely in pitch forward of its aerodynamic center, and pitch control is 
, achieved by a trimming surface, which is trimmed by trailing-edge flaps. The 
trimming surface is also free to pivot freely in pitch forward of its aerody- 
namic center but is connected to the wing by a rigidly attached boom. This 
report describes an analytical study of the longitudinal behavior of represen- 
tative small free-wing/ free- trimmer aircraft. It includes an assessment of the 
response to symmetric vertical turbulence, the nature of the characteristic 
longitudinal stick-fixed modes of motion, and the maximum trimmed lift 
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coefficients obtainable. Both forward- and aft-mounted trimming surfaces are 
considered. To avoid the lateral-directional instability problems found in 
previous free-wing studies, the left and right wing panels were constrained 
to symmetrical deflections only, roll control being achieved by the differen- 
tial operation of wing flaps. Body attitude adjustment was provided, as in the 
previous studies, by a separate horizontal stabilizer. 

Scope 

The research effort described in this report is limited to an analy- 
tical study of the stick-fixed longitudinal motion and gust alleviation charac- 
teristics of a representative light aircraft employing the f ree-wing/f ree- 
trimmer concept. The analysis was supported by the design of several concep- 
tual f ree-wing/free-trimmer aircraft to provide representative dimensional and 
mass parameters. 

Since both the free-wing and free-trimmer panels were constrained to 
symmetrical pitching motion, only the longitudinal behavior was analyzed. With 
no differential freedom permitted, the aircraft could be expected to behave 
similarly to conventional fixed-wing aircraft in the lateral-directional modes. 

Attention was confined to linear analyses of the responses to sym- 
metric vertical gusts and the basic longitudinal flying qualities as indicated 
by the characteristic roots associated with the stick-fixed longitudinal modes. 

The turbulence responses were evaluated by computing rms values of 
pertinent variables in response to continuous atmospheric turbulence, but no 
attempt was made for a direct evaluation of riding qualities as affected by 
human tolerance factors. 

The fundamental strategy of this investigation was to: 

(1) derive the equations of motion applicable to an aircraft using 
the f ree-wing/free-trimmer concept ; 

(2) create a family of conceptual aircraft designs to provide realis 
tic parameters for the mathematical models; and 

(3) perform a linear analysis of the gust alleviation character- 
istics, handling qualities, and other features peculiar to the overall ‘ concept . 
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mean aerodynamic chord, meters (feet) 

maximum section lift coefficient 

free-trimmer lift coefficient 

wing lift coefficient 

maximum wing lift coefficient 

moment coefficient about aerodynamic center 

moment coefficient about hinge axis 

components of force transmitted from free trimmer to boom, newtons 
(pounds) 

components of force transmitted from wing to fuselage, newtons 
(pounds) 

hinge margin, distance from aerodynamic center forward to hinge 
axis, expressed as fraction of chord length 

wing reference area, meters^ (feet^) 

free-trimmer reference area, meters^ (feet^) 

increment in airspeed divided by equilibrium airspeed 

distance of free-trimmer hinge forward of wing hinge, meters (feet) 

angle of attack of free-trimmer surface, radians 

angle of attack of fuselage assembly, radians 

angle of attack of wing, radians 

angular displacement of free trimmer with respect to wing, radians 
wing flap deflection, degrees 

angular displacement of wing with respect to fuselage, radians 
pitch angle of fuselage with respect to horizon, radians 
rms normal load factor response, g ! s 
Laplace operator, 1/sec 


PROCEDURE 


Mathematical Models 

Although the original plan was to modify the computer programs used 
in the previous free-wing analyses, it soon became apparent that the additional 
complexity introduced by the free-trimmer dynamics dictated the development of 
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new mathematical models. Furthermore, whereas the linear models used in the i 
previous analyses were derived directly from the complete nonlinear equations 

i 

by limiting the motion to small disturbances from straight and level equilib- 
rium, the current task required an unorthodox form of the linear model. | 

The development of the inertial force and moment terms for the com- 

i 

plete nonlinear equations with eight degrees of freedom is given in appendix A , 
to provide documentation for possible future studies. The linear development, 
which is confined to longitudinal motion, is described separately in appendix B. 

The longitudinal system has five degrees of freedom. Three variables 
are required to define the spatial position and orientation of the fuselage in 
the vertical plane, and two additional variables are required to define the 
angular displacements of the wing and trimmer surfaces. With unsteady aero- 
dynamic effects included for both lifting surfaces, the linear longitudinal 
system is of tenth order*. 

The linearized set of equations describing the longitudinal motion 
of the aircraft in response to vertical gust velocities and control tab dis- 
placement is given by the following equation: 

V 

[A] X = [B] 6 1 + [C] ^ (1) 

u o 

Where [A] is a 13 by 13 matrix of the coefficients of the homogeneous equations, 

and [B] and [C] are column matrices. 

The state variable X has 13 components: a c > , 6^, F x , F z , a , 

c c-*b c->b 


C L 5 6 p’ F x ’ F z ’ 6 ’ u ’ and a f‘ 
w F w->f w-*f 

It is the number of components of the state vector which constitutes 
the primary departure from convention. Since the system has five independent 
degrees of freedom, it would have been possible to combine most of the equa- 
tions to arrive at a set using only five independent variables. These would 
logically have been u, a^, 0, 6 , and 5^ . Unfortunately, the elimination of 
the eight remaining variables involves considerable algebraic manipulation 
which would have obscured the physical significance and origins of the indivi- 
dual terms in the final equations. 


* The linear approximations of reference 5 were used to model the unsteady lift 
phenomena, as discussed in appendix B. 
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Although the retention of 13 variables in the equations of motion 
seems ponderous, and incurs certain penalties in computation time, it has the 
overriding advantage of greatly facilitating the verification of each of the 
terms in the individual equations — an advantage of some importance, since no 
prior analysis of such a system was in existence to provide a standard for 
checking the correctness of the model. Since time did not permit the restruc- 
turing of the model in condensed, more conventional form, the 13 equations were 
retained as they were derived. 

Conceptual Aircraft Designs 

In this study, attention was confined to families of airplanes based 
upon two light aircraft configurations in order to assess some benefits, if 
any, accruing to light aircraft of markedly different missions. The first was 
a two seat light aircraft having tricycle gear, pusher propeller, and a high 
free wing with a free trimmer placed two chord lengths ahead of the wing hinge 
line and fixed to the wing by two rigid booms. The second aircraft was a 
single seat agricultural airplane having conventional gear and a tractor pro- 
peller. The free wing was mounted lower on the fuselage and the free trimmers 
were mounted aft of the wing on two rigid booms. 

For purposes of direct comparison and parametric analysis of the 
penalty and benefit of free trimmer control volume, both forward and aft of 
the wing, the first configuration was fitted with a forward free trimming 
surface, referred to in this work as a forward trimmer, two separate aft free 
trimming surfaces referred to as aft trimmers, a pure free wing, and a compara- 
ble fixed wing. 

Design summary. - Following accepted procedures for the conceptual 
design of light aircraft, the two types of aircraft were designed in enough 
detail to allow the calculation of weights and moments of inertia. Steps in 
this process common to both were the identification of aircraft mission, the 
calculation of the flying surface sizes necessary to fulfill their missions, a 
first weight estimate based on mission requirements, an estimate of wing maxi- 
mum trimmed lift, tail sizing, development of a cockpit design that conformed 
to MIL-STD-1333 , first layout, revision of the weight estimate, iterations to 
reduce weight, and the calculation of moments of inertia. Common features to 
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both aircraft are an aspect ratio 6 rectangular wing with a NACA 23012 section, 
wing loading comparable to existing aircraft of the same class, a total free 
trimmer area that is one-sixth the wing area, a horizontal tail sized strictly 
for fuselage pitch attitude control, wing panels and trimmer panels that are 
constrained to symmetrical deflections only, and 19 percent chord hinge lines 
for all free surfaces. 

Previous work dictated that the f ree-wing/boom/trimmer combination 
balance about the hinge point of the wing. This is an important design con- 
straint which was maintained in the conceptual configurations by carrying all 
disposable items in the fuselage. Variations of the center of gravity of the 
trimmer, wing/boom/ trimmer , and fuselage were permitted for parametric analyses 
during the investigation. 

First concept - a two place light aircraft. - This is a two place 
aircraft with the free trimmer in a forward position ahead of its free wing. 

The basic mission for sizing is fair weather training. Payload is similar to 
existing aircraft in this mission category and includes two crew and enough 
fuel for a 3 hour flight. Fixed wing gross weight is in the 680 kilogram 
(1500 pound) class. Other features are two place side by side seating ahead 
of the wing, a shoulder wing, typical general aviation wing loading, an 
86 kilowatt (115 BHP) engine in a pusher arrangement, 20 percent chord flaps on 
the free wing and free trimmer, and a two wing-chord-length arm on the free 
trimmer. 

The family of configurations is listed in table I and figures 2(a) to 

2(c) are multiview drawings of each of the parametric variations described in 
table I below: 


TABLE I. FAMILY OF TWO-PLACE LIGHT AIRCRAFT 


Configuration 

Figure 

Description 

1A 

2(a) 

Forward free surface 

IB 

2(b) 

Pure free wing 

1C 

2(c) 

Aft free surface 

ID 

(Not shown) 

Fixed-wing equivalent to 
configuration IB 
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FIGURE 2(b). TWO SEAT LIGHT AIRCRAFT (CONFIGURATION IB- 

PURE FREE WING) . FIXED WING CONFIGURATION ID 
EXTERNALLY IS SIMILAR TO THIS LESS EXTERNAL 
BALANCES ON WING. 






Table II is a summary of geometric data, and table III is a weight ^ 
summary for each of the configurations in table I. No attempt was made to 
minimize the counterbalancing weight penalty due to static balance require- 
ments to provide a direct comparison of' benefits and penalties from one | 

configuration to another. j 

These configurations are strictly conceptual and address only the 

i 

most obvious practical design constraints. The propeller and engine place- 
ment, for instance, may not be practical in a real aircraft because of lack of 
propeller ground clearance. 

Second concept - a single place agricultural aircraft. - This air- 

! 

craft was sized to be comparable to existing American agricultural aircraft in 
the 1700 kilogram (3800 pound) weight range with about a 0.5 cubic meter 
(17 cubic foot) hopper capacity. Features of this aircraft design are a single j 

place air conditioned and pressurized cockpit surrounded by 40g crash-resis- 
tant shell, a 0.45 cubic meter (16 cubic foot) hopper (approximately 590 

t 

kilogram (1300 lb) payload) , a 227 kilowatt (305 BHP) engine with a constant speedy 
propeller, a fully cantilevered wing, spray bars carried in external struts 

under each wing and hinged on the wing main spar, conventional gear, and a mid- 
wing to eliminate the need for geometric dihedral. 

The family of parametric configurations examined is shown in 
figures 3(a) and 3(c) and their designations are explained in table IV. 

Geometric data are presented in table V, and weight data are shown in table VI. 
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TABLE II. SUMMARY OF GEOMETRIC DATA FOR TWO SEAT 
AIRCRAFT VARIATIONS 


Item 

1A 

Fuselage - 


Length, m (ft) 

6.44 (21.13) 

Maximum width, m (ft) 

1.19 (3.9) 

Maximum height, tn (ft) ! 

1.31 (4.3) 

Surface area, m 2 (ft 2 ) j 

18.247 (196.4) 

Wing - j 


Span, m (ft) 

7.77 (25.5) 

Reference area. 


m 2 (ft 2 ) i 

10.09 (108.6) 

Exposed area, m 2 (ft 2 ) 

8.54 (91.9) 

Wetted area, m 2 (ft 2 ) 

17.38 (187.1) 

Aspect ratio i 

6.00 

Taper ratio ! 

1.00 

Thickness ratio \ 

0.12 

Mean geometric 


chord, m (ft) 

1.30 (4.25) 

Airfoil ! 

NACA 23012 

Hinge point (percent 


chord) | 

19.0 

Free stabilizer (each 


surface if more than 


one) - j 


Span, m (ft) j 

3.17 (10.4) 

Reference area , 


m 2 (ft 2 ) 

1.68 (18.1) 

Exposed area, :n 2 (ft 2 ) 

1.68 (18.1) 

Wetted area, U' 2 (ft^) j 

3.42 (36.8) 

Aspect ratio \ 

6.0 

Taper ratio j 

1.0 

Thickness rati 5 

0.12 

Mean geometric i 

i 

chord, m (ft) 

0.52 (1.7) 

Airfoil 

j NACA 23012 

Hinge point (percent 

i 

chord) 

119.00 

Control arm (hinge 

i 

i 

point to hinge 

i 

point) , m (ft) 

1-2.53 (-8.3) 

Control volume 


(hinge point to 

i 

hinge point) 

j-0.33 

Horizontal tail - 


Span, m (ft) 

2.44 (8.0) 

Reference area. 


m 2 (ft 2 ) 

2.14 (23.0) 

Exposed area, m 2 (ft 2 ) 

1.49 (16.0) 


Configuration 

IB . 1C ID 


6.44 

(21.13) 

1.19 

(3.9) 

1.31 

(4.3) 

18.247 

(196.4) 

7.77' 

(25.5) 

10.09 

(108.6) 

8.54 

(91.9) 

17.38 

(187.1) 

6.00 


1.00 


0.12 


1.30 

(4.25) 

NACA 

23012 


19.0 


2.44 (8.0) 

2.14 (23.0) 
1.49 (16.0) 


6.44 

(21.13) 

1.19 

(3.9) 

1.31 

(4.3) 

18.247 

(196.4) 

7.77 

(25.5) 

10.09 

(108.6) 

8.54 

(91.9) 

17.38 

(187.1) 

6.00 


1.00 


0.12 


1.30 

(4.25) 

NACA 

23012 

19.0 



1.14 

(3.8) 

0.871 

(9.4) 

0.871 

(9.4) 

1.77 

(19.1) 

1.50 


1.00 


0.12 


0.76 

(2.5) 

NACA 

23012 

19.00 


1.31 

(4.3) 

0.18 


2.44 

(8.0) 

2.14 

(23.0) 

1.49 

(16.0) 


6.44 

(21.13) 

1.19 

(3.9) 

1.31 

(4.3) 

18.247 

(196.4) 

7.77 

(25.5) 

10.09 

(108.6) 

8.54 

(91.9) 

17.38 

(187.1) 

6.00 


1.00 


0.12 


1.30 

(4.25) 

NACA 

23012 


None 


2.44 (8.0) 

2.14 (23.0) 
1.49 (16.0) 
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TABLE II. (Continued) 





Configuration 



Item 

1A 


IB 

1C 


ID 

Wetted area, m2 (ft2) 

3.03 (32.6) 

3.03 

(32.6) 

3.03 (32.6) 

3.03 

(32.6) 

Aspect ratio 

2.78 

2.78 


2.78 

2.78 


Taper ratio 

0.50 

0.50 


0.50 

0.50 


Thickness ratio 
Mean geometric 

0.09 

0.09 


0.09 

0.09 


chord, m (ft) 

0.951 (3.12) 

0.951 (3.12) 

0.951 (3.12) 

0.951 (3.12) 

Airfoil 

Control arm (hinge 

NACA 0009 

NACA 

0009 

NACA 0009 

NACA 

0009 

point to one-quarter 
chord), m (ft) 

3.05 (10.0) 

3.05 

(10.0) 

3.05 (10.0) 

3.05 

(10.0) 

Control volume (hinge 







point to one-quarter 







chord) 

0.50 

0.50 


0.50 

0.50 


Vertical tail - 







Span, m (ft) 
Reference area 

1.02 (3.3) 

1.02 

(3.3) 

1.02 (3.3) 

1.02 

(3.3) 

m2 (ft2) 

2.02 (21.7) 

2.02 

(21.7) 

2.02 (21.7) 

2.02 

(21.7) 

Exposed area m2 (ft2) 

2.02 (21.7) 

2.02 

(21.7) 

2.02 (21.7) 

2.02 

(21.7) 

Wetted area m2 (ft2) 

4.097 (44.1) 

4.097 (44.1) 

4.097 (44.1) 

4.097 (44.1) 

Aspect ratio 

0.51 

0.51 


0.51 

0.51 


Taper ratio 

0.50 

0.50 


0.50 

0.50 


Thickness ratio 
Mean geometric chord. 

0.09 

0.39 


0.09 

0.09 


m (ft) 

1.04 (3.42) 

1.04 

(3.42) 

1.04 (3.42) 

1.04 

(3.42) 

Airfoil 

Control £ rm (hinge 

NACA 0009 

NACA 

0009 

NACA 0009 

NACA 

0009 

point to one- 
quarter chord) , m 
(ft) 

2.92 (9.6) 

2.92 

(9.6) 

2.92 (9.6) 

2.92 

(9.6) 

Control volume (hinge 







point to one- 







quarte. chord) 

Pitch moment of inertia ■ 
Forward trimmer about 
trimmer hinge line, 
kg-m 2 (slug-f t2) 

Aft trimmer about 

0.08 

1.2 (1.637) 

0.08 


0.08 

0.08 




36.0 (48.839) 



trimmer hinge line, 
kg-m 2 (slug-ft 2 ) 


129.7 (175.894) 

34.1 

(46.283) 



Wing/boom/balances 
about wing hinge 
line, kg-m2 (slug 
ft 2 ) 

56.8 (77.027) 

77Q i | 

(1056.35) 

Total aircraft about 
aircraft center of 
gravity, kg-m 2 
(slug-ft 2 ) 









/J3.1 1 
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TABLE III. 


Item 

1A 

Total Structure - 

353 

(779) 

Wing group including 



balances 

: 148 

(326) 

Tail group: 



Horizontal 

10 

(21) 

Vertical 

5 

(ID 

Body group 

88 

(193) 

Alighting gear group: 



Main 

24 

(54) 

Auxiliary 

5 

(10) 

Surface controls 

21 

(47) 


Engine section 

Booms 15 (33) 


Trimmer including 
balances 


Fixed Equipment - 
Instruments 
Electrical group 
Furnishings 
Total Weight Empty - 
Crew 
Fuel: 

Internal 

Trapped 

Oil: 


38 (84) 


| Propulsion Group - 

156 

(345) 

Engine (as installed) 

: 111 

(244) 

Accessories 

2 

(5) 

Air induction system 

5 

(10) 

Exhaust system 

1 

(2) 

Cooling and drain 



provisions 

1 

(3) 

Lubricating system 

1 

(2) 

Fuel system 

7 

(15) 

Engine controls 

1 

(2) 

Starting system 

11 

(24) 

Propeller installation 

12 

(26) 

Drive shaft 

5 

(12) 


54 (120) 
3 (6) 


Engine 

Trapped 

Total Useful Load 
Takeoff Gross Weight 
Flight Design Gross Weight 
Landing Design Gross Weight 


3 (7) 
0.5 (1) 
215 (474 


783 (1727) 
761 (1679) 


740 (1631) 


SUMMARY 
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FIGURE 3(b). 


SINGLE SEAT AGRICULTURAL AIRCRAFT (CONFIGURATION 2B — PURE 
FREE WING) . FIXED WING CONFIGURATION 2D IS EXTERNALLY 
SIMILAR TO THIS LESS EXTERNAL BALANCES ON WING. 




TABLE IV. FAMILY OF AGRICULTURAL AIRCRAFT 


Configuration 

Figure 

Description 

2A 

3(a) 

Forward free surface 

2B 

3(b) 

Pure free wing 

2C 

3(c) 

Aft free surface 

2D 

Not shown 

Fixed-wing equivalent to 
Configuration 2B 


TABLE V. SUMMARY OF GEOMETRIC DATA FOR SINGLE SEAT 
AGRICULTURAL AIRCRAFT VARIATIONS 


Item 

Configuration 

2A 2B 

2C 

2D 

Fuselage - 

1 

f 


Length, m (ft) 

7.5 (24.50) ! 7.5 (24.50) 

7.5 (24.50) ! 

7.5 (24.50) 

Maximum width, m (ft) 

j 1.0 (3.25) | 1.0 (3.25) i 

1.0 (3.25) ! 

1.0 (3.25) 

Maximum height, m (ft) 

| 1.6 (5.25) | 1.6 (5.25) j 

1.6 (5.25) 

1.6 (5.25) 

Surface area, (ft 2 ) 

j20.1 (215.92) ; 20.1 (215 . 92); 20. 1 (215.92) 

20.1 (215.92; 

Wing - 

! 



Span, m (ft) 

110.5 (34.6) 10.5 (34.6) 

10.5 (34.6) 

10.5 (34.6) 

Reference area. 

S i 



m 2 (ft 2 ) 

19.3 (207.5) 19.3 (207.5) 

19.3 (207.5) 

19.3 (207.5) 

Exposed area, m 2 (ft 2 ) 

17.5 (188.8) 17.5 (188.8) 

17.5 (188.8) 

17.5 (188.8) 

Wetted area, m 2 (ft 2 ) 

|35.7 (384.1) 35.7 (384.1) 

35.7 (384.1) 

35.7 (384.1) 

Aspect ratio 

| 6.0 6.0 

6.0 

6.0 

Taper ratio 

I 1.0 ; 1.0 

1.0 

1.0 

Thickness ratio 

! 0.12 ; 0.12 

! 0.12 

0.12 

Mean geometric chord. 




m (ft) 

1.8 (5.76) 1 1.8 (5.76) 

1.8 (5.76) 

1.8 (5.76) 

Airfoil 

NACA 23012 ; NACA 23012 

NACA 23012 

NACA 23012 

Hinge point (percent 

| 



chord) 

19.0 19.0 

19.0 

19.0 

Free stabilizer (each 

t 

1 



surface if more than 

; [ 
i 



one) - 




Span, m (ft) 

2.2 (7.1) | 

1.5 (5.0) 


Reference area, 

I 



m^ (ft2) 

1.6 (16.8) 

1.6 (16.7) 


Exposed area, m2 (ft 2 ) 

1.4 (15.2) j 

1.6 (16.7) 


Wetted area, m 2 (ft 2 ) 

; 2.9 (31.0) j 

3.1 (33.9) 


Aspect ratio 

3.0 ! 

1.5 

j 

Taper ratio 

1.0 ! 

1.0 

! 

Thickness ratio 

! 0.12 | 

0.12 

! " 

Mean geometric chord, 

! i 

i 

| 

m (ft) 

0.7 (2.4) ! 

1 1.0 (3.3) 

| 
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TABLE V. (Continued) 






Configuration 



Item 


2A 


2B 


2C 


2D 

Airfoil 

NACA 

23012 

j „ 

« 

! NACA 

23012 




Hinge point (percent 









chord) 

19.0 


— 


19.000 

— 


Control arm (hinge 









point to hinge 









point), m (ft) 

-1.7 

(-5.5) 

— 


: 1.8 

(5.8) 

— 


Control volume (hinge 









point to hinge 









point) 

-0.15 




0.16 




Horizontal tail - 









Span, m (ft) 

3.1 

(10.3) 

3.1 

(10.3) 

3.1 

(10.3) 

3.1 

(10.3) 

Reference area. 









m^ (ft^) 

3.2 

(34.6) 

3.2 

(34.6) 

3.2 

(34.6) 

3.2 

(34.6) 

Exposed area, m^ (ft^) 

2.9 

(31.3) 

2.9 

(31.3) 

2.9 

(31.3) 

2.9 

(31.3) 

Wetted area, m^ (ft^) 

5.9 

(63.7) 

5.9 

(63.7) 

5.9 

(63.7) 

5.9 

(63.7) 

Aspect ratio 

3.0 


3.0 


3.0 


3.0 


Taper ratio 

0.5 


0.5 


0.5 


0.5 


Thickness ratio 

0.09 


0.09 


0.09 


0.09 


Mean geometric chord 

1.1 

(3.5) 

1.1 

(3.5) 

1.1 

(3.5) 

1.1 

(3.5) 

Airfoil 

NACA 

0009 

NACA 

0009 

NACA 

0009 

NACA 

0009 

Control arm (hinge 









point to one-quarter 









chord) , m (ft) 

4.2 

(13.92) 

4.2 

(13.92) 

4.2 

(13.92) 

4.2 

(13.92) 

Control volume (hinge 









point to one-quarter 









chord) 

0.4 


0.4 


0.4 




Vertical tail - 









Span, m (ft) 

1.1 

(3.5) 

1.1 

(3.5) 

1.1 

(3.5) 

1.1 

(3.5) 

Reference area, 









m^ (ft^) 

2.1 

(22.5) 

2.1 

(22.5) 

2.1 

(22.5) 

2.1 

(22.5) 

Exposed area, m^ (ft^) 

2.1 

(22.5) 

2.1 

(22.5) 

2.1 

(22.5) 

2.1 

(22.5) 

Wetted area, m^ (ft^) 

4.2 

(45.7) 

4.2 

(45.7) 

4.2 

(45.7) 

4.2 

(45.7) 

Aspect ratio 

0.5 


0.5 


0.5 


0.5 


Taper ratio 

0.5 


0.5 


0.5 


0.5 


Thickness ratio 

0.09 


0.09 


0.09 


0.09 


Mean geometric 









chord, m (ft) 

1 1.1 

(3.5) 

1.1 

(3.5) 

1.1 

(3.5) 

1.1 

(3.5) 

Airfoil 

| NACA 

. 0009 

NACA 

i 

. 0009 

J 

NACA 

. 0009 

NACA 

L 

. 0009 
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TABLE VI. WEIGHT SUMMARY 


Total Structure - 

Wing group including 

balances 

Tail group: 

Horizontal 
Vertical 
Body group 

Alighting gear group: 

Main 

Auxiliary 
Surface controls 
Engine section 
Booms 

Trimmer including 
balances 

PrcDulsion Group - 

Engine (as installed) 
Accessories 
Air induction system 
Exhaust system 
Cooling & drain provisions 
Lubricating system 
Fuel system 
Engine controls 
Jtarting system 
Propeller installation 
Fiy id Equipment - 
Instruments 
Electrical group 
Electronics group 
F irnishings 
Air conditioning 
equipment 

Auxiliary gear (spray 
tank and nozzles) 

Total Weight Empty - 
Crew 
Fuel: 

Internal 

Trapped 

Oils 

Engine 

Trapped 

Equipment (AG spray 
or dust) 

Total Useful Load 
Takeoff Gross Weight 
Flight Design Gross Weight 
Landing Design Gross Weight 


Configuration 

2A 2B 2C 2P 

537 (1184) I 586 (1291) 654 (1441) 465 (1026) 


25 (56) 

7 (16) 
121 (266) 

64 (141) 
5 (11) 
24 (54) 
10 ( 21 ) 
( 68 ) 

(90) 

268 (591) 


25 (56) 

7 (16) 
121 (266) 

64 (141) 
5 (11) 
24 (54) 
10 ( 21 ) 


268 (591) 


154 (340) 


154 (340) 


25 (56) 

7 (16) 
121 (266) 

64 (141) 

5 (11) 

24 (54) 

10 ( 21 ) 
( 112 ) 

(90) 

268 (591) 

205 (452) 

3 (7) 

1 ( 2 ) 

1 ( 2 ) 

1 ( 2 ) 

2 (5) 

10 ( 22 ) 

1 ( 2 ) 

11 (24) 

33 (73) 
(154 (340)) 

18 (40) 
11 (25) 
16 (35) 
23 (51) 

34 (76) 


25 (56) 

7 (16) 
121 (266) 

64 (141) 
5 (11) 
24 (54) 
10 ( 21 ) 


268 (591) 


154 (340) 


j I 51 (113) 

959 (2115) : 1008 (2222) 1076 (2372) 

i i 77 (170) 


888 (1957) 


106 (234) 
3 (6) 

6 (14) 

1 ( 2 ) 


590 (1300) 
783 (1726) - 783 (1726) ; 783 (1726) 
1742 (3841) i 1790 (3948) ; 1859 (4098) 
11699 (3747) ! 1748 (3854) (1816 (4004) 
i 1657 (3654) i 1706 (3761) 11774 (3911) 


783 (1726) 
1670 (3683) 
i 1628 (3589) 
11585 (3496) 


23 


Handling Qualities 


Handling qualities criteria were based upon the requirements of the 
revised military handling qualities specification, MIL-F-8785B (ASG) , and were 
confined to an examination of the damping of the phugoid and longitudinal short- 
period modes for stick-fixed motion. Compliance with the standards was deter- 
mined by examining the characteristic roots of the system of equations. It is 
recognized that the applicability of the conventional handling qualities cri- 
teria has not been established for free-wing aircraft. 

When represented as a polynomial in the operator X, the determinant 
of the matrix of coefficients of equation (1) becomes a tenth degree expression 
with eleven polynomial coefficients. To compute these coefficients, eleven 
arbitrary real values of the operator were selected, and the corresponding 
values of the determinant were computed. Then a set of eleven simultaneous 
equations was formed and solved for the eleven coefficients. The roots of this 
polynomial were then obtained by a standard computer subroutine to arrive at 
the characteristic roots of the system. The 10 roots were typically divided 
into four complex pairs and two heavily damped real roots. 

Responses to Atmospheric Turbulence 
The responses to atmospheric turbulence were computed by using the 
techniques outlined in appendix C. A power spectral density technique was used 
to obtain the statistical responses to continuous turbulence, and a discrete 
gust technique was used to examine the details of an encounter with an isolated 
gust disturbance. Only vertical gust components were considered. 

DISCUSSION OF RESULTS 

Longitudinal Characteristic Modes 
In this section, the characteristic modes of a free-wing/free- 
trimmer aircraft are discussed by the use of specific numerical examples, and 
comparisons are made with the natural motions of conventional fixed-wing air- 
craft and the pure free-wing aircraft analyzed in references 1 ana 2. 

It is well known that the longitudinal motion of a conventional rigid 
aircraft with controls fixed is adequately described by a set of equations 
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' yielding four characteristic roots. These four roots are typically divided 
I into two complex pairs, one pair defining the long-period phugoid motion, and 

! the other pair representing the longitudinal short-period mode, 

j For a pure free-wing aircraft with trim provided by trailing-edge 

r control surfaces, the previous analyses have shown that an additional oscilla- 
tory mode appears. The new oscillatory mode, which was termed the symmetric 
wing panel mode in reference 1, was easy to distinguish from the usual short- 
period mode because of its higher frequency of oscillation (a consequence of 
the low pitching inertia of the wing panels) , and it was the dominant mode 
with regard to normal load factor responses to control inputs and vertical gust 
1 excitation. The roots associated with the conventional short-period mode, on 

! 

^ the other hand, described a motion almost totally confined to the pitching of 
the fuselage assembly. From an analysis standpoint, the clear physical and 
1 numerical separation of these modes permitted a relatively straightforward 

assessment of the significance of parametric changes in the motions of the 
aircraft. 

In the current analysis, the addition of a separate free trimming 
surface creates yet another oscillatory mode that is related to the pitching 
motion of the trimmer surface. Although the new mode is easily identified by 
its high natural frequency and has a readily understood physical significance, 
the addition of the separate surface complicates the analysis because of the 
much greater coupling which it introduces between the motions of the fuselage 
and wing/trimmer assemblies. Not only does the increased pitching inertia of 
the wing/trimmer place the natural frequencies of the wing and fuselage into 
closer proximity; but in addition, the transient aerodynamic trimmer forces 
caused by the vertical motion of the entire aircraft create significant pitch- 
ing moments about the wing hinge axis. As a consequence, in the terminology of 
reference 1, the symmetric wing panel mode and the short-period mode cannot, in 
general, be easily distinguished with regard to their physical manifestation 
in the behavior of the aircraft. For this reason, the two modes are simply 
termed mode A and mode B in the discussion which follows. 

The set of linearized equations derived in appendix B defines a tenth 
order dynamic system for controls-f ixed motion. Typically, the 10 character- 
istic roots include four complex pairs denoting oscillatory modes and two 
heavily damped aperiodic modes related to unsteady aerodynamic forces. 
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To provide an illustration of modal behavior, table VII lists the 
numerical values of the characteristic roots for two configurations which 
differ primarily in the location of the trimming surface. The first (forward 
trimmer, configuration 1A, fig. 2(a)) has a trimming surface located two wing 
chord lengths forward of the wing hinge axis. The second configuration has 
the same trimmer area and moment arm, but the trimmer is aft of the wing hinge. 
Both configurations have a mass balanced fuselage, wing assembly, and trimmer 
about their respective hinge axes. Calculations were made for a normal cruise 
condition of 241.4 kilometers per hour (150 miles per hour) at an altitude of 
1829 meters (6000 feet) . 


TABLE VII. 


COMPARISON OF CHARACTERISTIC ROOTS FOR FORWARD 
AND AFT TRIMMING SURFACES 


=2.0, CRUISE FLIGHT 


c 


Mode 

Forward trimmer 

Aft trimmer 

Phugoid 

-0.0509 + j 0.321 

-0.0676 + j 0.313 

Mode A 

-0.171 + j 2.42 

-1.890 + j 2.35 

Mode B 

-2.38 + j 5.12 

-0.119 + j 4.42 

Trimmer mode 

-6.60 + j 17.2 

-7.57 + j 17.7 

Aperiodic 

-30.4, -75.3 

-30.5, -75.1 


Figure 4 i s a comparison of the responses to step displacements, of 
equal magnitude, of trailing-edge control flaps on the trimmer surfaces, and 
illustrates the inadvisability of attempting to assign a simple physical sig- 
nificance to either mode A or mode B. 

Even though modes A and B are the respective generic equivalents of 
the short-period and symmetric wing panel modes identified in references 1 and 
2, they are clearly not uniquely related to either the fuselage or to wing mo- 
tion alone. One of these modes is much more heavily damped than the other, but 
the relative damping is reversed between the two configurations. Furthermore, 
the mode with the lighter damping ratio appears to dominate the response. The 
forward trimmer response displays the relatively long period of mode A, while 
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the aft trimmer response is governed by the shorter period of mode B; the 
dominant mode in either case is the one with the lighter damping. 

It should be noted that these cases are shown only for illustration, 
and it is recognized that neither represents acceptable response characteristics 
from a handling qualities standpoint. The shortcomings of configurations such 
as these, with large trimmer moment arms, will be discussed in the following 
sections . 


Configuration Options and Limitations 

The most fundamental design feature of the free-wing/free-trimmer 
concept is the size and location of the trimming surface relative to the wing. 

A key portion of this investigation was therefore an assessment of the gross 
effects of wing/trimmer geometry. 

For this evaluation, the pure free-wing aircraft version of the light 
aircraft of figure 2(b) (configuration IB) was used as a baseline. The fuse- 
lage and basic wing were unchanged, and the trimming surface and its supporting 
structure, as well as the necessary counterweights, were added to the wing. 

For this portion of the analysis, mass balance was maintained about the hinge 
axes of the trimmer surface and the wing /boom/ trimmer assembly. In this way, 
the general effects of trimmer geometry on wing assembly moment of inertia were 
implicitly included in the parametric changes, and the relative weight pen- 
alties for the various configurations could be evaluated. 

In assessing the merits of each configuration, four criteria were 
used: gust alleviation, handling qualities, maximum obtainable trimmed lift 
coefficient, and relative weight penalty incurred by mass balancing. 

The gust alleviation characteristics of each configuration were 
evaluated by comparing its vertical gust response to the vertical gust response 
of the baseline pure free-wing aircraft and to a geometrically similar fixed 
wing airplane. The comparison criteria were the rms normal load factor re- 
sponse to continuous vertical turbulence and the peak load factor values com- 
puted for an encounter with a discrete 1-cosine gust with a wavelength of 
25 wing chords and , ak magnitude of 3.048 meters per second (10 feet per 
second) . 

For the continuous turbulence analysis, the velocity degree of free^ 
dom was eliminated from the equations of motion to avoid large contributions 
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at very low frequencies near the phugoid mode, since these are easily controlled 
by pilot action. The Dryden turbulence spectrum was used with a scale length 
of 533 meters (1750 feet) for consistency with the previous work of references 
1 and 2. Details of the procedure for computing turbulence responses are con- 
tained in appendix C. 

Although the continuous turbulence analysis is probably the more 
meaningful of the two, the discrete gust calculations, in time history form, 
provided valuable insight into the detailed behavior of the aircraft during 
gust encounters. The 25 chord length gust was used to conform to the stan- 
dards of the Federal Air Regulations (FAR 23 . 333(2) (2) (i) ) , which specify this 
type of discrete gust for analytical proof of compliance with the structural 
gust load criterion. A comparison of the peak discrete gust loads with the 
loads of the equivalent fixed-wing aircraft therefore yields a rapid, first 
order estimate of the relative structural gust loads which each aircraft would 
be required to withstand for certification purposes. 

With regard to handling qualities, the primary criterion was the 
damping ratio of the more poorly damped of the two short-period modes. The 
military handling qualities criterion, MIL F-8785B(ASG) , requires the damping 
ratio of the short-period motion to lie between 0.35 and 1.30 for level 1 per- 
formance.* For level 3 flying qualities, a damping ratio of 0.15 is permitted. 

Forward trimmer configurations. - Since the trimming force with for- 
ward trimming surfaces is in an upward direction, thereby augmenting wing lift, 
the forward trimmer configurations were initially expected to offer the most 
promise. The investigation began, therefore, with an analysis of Configura- 
tion 1A (fig. 2(a)), with a forward trimmer moment arm of two wing chord 
lengths and an area ratio of 1/6. 

*MIL-F-8785B (ASG) defines three levels of acceptability: 

Level 1: Flying qualities clearly adequate. 

Level 2: Flying qualities adequate to accomplish the mission... 

but some increase in pilot work load or degradation 
of mission ef fectiveness exists. 

Level 3: Flying qualities such that the airplane can be con- 

trolled safely but pilot workload is excessive or 
mission effectiveness is inadequate, or both. 
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The response to continuous turbulence of unit intensity was computed 
for the cruise condition of 241.4 kilometers per hour (150 miles per hour) at an 
altitude of 1829 meters (6000 feet). The response proved to be disappointing - 
higher than the equivalent fixed-wing aircraft. Table VIII compares the rms 
responses to continuous turbulence as well as the peak load factor responses to 
discrete gusts for the f ree-wing/f orward free trimmer of Configuration 1A, the 
pure free-wing version of Configuration IB (fig. 2(b)), and the fixed wing 
equivalent, Configuration ID. The input magnitude used for the Dryden model 
of continuous turbulence is not related to the discrete gust input magnitude. 
Instead, a value of 1 foot per second, or unity, was selected as the rms value 
for reasons of simplicity. As a consequence, the rms load factors shown are 
appreciably less than would actually be measured in turbulence. 


TABLE VIII. COMPARISON OF TURBULENCE RESPONSES 


Configuration 

Rms load factor 
per unit turbulence 
intensity 

Peak load factor 
response to 
discrete gust 

Free wing/ 



free forward trimmer 



(conf. 1A) 

0.0273 

0.730 

Free wing (conf. IB) 

0.00989 

0.348 

Fixed wing (conf. ID) 

0.0191 

0.675 


The poor gust response of the forward free-trimmer configuration can 
be explained by an examination of the response histories in figure 5. For 
these encounters, all three of the aircraft of Table VIII were subjected to the 
standard 25 chord-length 1 - cosine gust velocities shown at the bottom of the 
figure. In the pure free-wing version, the wing begins to deflect almost 
immediately after encountering the gust, thereby limiting the load factor 
response to the relatively low value shown in the upper trace. In contrast, the 
forward trimmer surface experiences an initial upward transient lift force 
before its downward pitching motion can relieve the load. The transient for- 
ward trimmer force, though lasting less than one-fifth of a second, imparts a 
pitch-up impuse to the wing assembly which not only delays the alleviating 
motion of the wing, but actually produces a small wing displacement in an 
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adverse direction. In the case shown, the wing does not begin to pitch down- 
ward until the gust has already passed its peak value. 

Aside from the poor gust response characteristics, Configuration 1A 
is unsatisfactory from a handling qualities standpoint because of the low damp- 
ing of the dominant short-period oscillation. As noted in figure 4, the 
dominant short-period mode is the oscillation identified as mode A in table VII, 
and the damping ratio of this mode is only 0.0705, which is well below the 
lowest standard (level 3) of MIL-F-8785B (ASG) . 

The poor gust response of Configuration 1A is related to the trimmer 
volume (area times moment arm) , while the low damping of the oscillatory mode 
may be a consequence of the relatively large rotational inertia of the wing/ 
boom/forward trimmer system as compared to the small pitch damping provided by 
the free forward trimmer at the modal frequency. 

To improve the behavior of the forward trimmer configurations, two 
changes were explored: the use of smaller moment arms and the relocation of the 
hinge axis on the forward trimmer surface. 

With regard to the hinge axis, the nominal hinge margin for the 
forward trimmer surface of Configuration 1A is 6 percent; that is, the hinge 
axis is located 6 percent of the trimmer chord forward of the aerodynamic 
center. From earlier free-wing analyses it is known that the greater the hinge 
margin, the more rapid the alleviation of out-of-trim conditions. It might be 
expected that the unfavorable pitch impulse imparted by the trimmer, as in the 
gust encounter in figure 5, would be reduced by quickening the response of the 
trimmer surface. On the other hand, the larger the hinge margin the greater 
the control tab size and deflection required for trim, and the smaller the max- 
imum trimmed lift coefficient of the free surface. In light of these factors, 
all subsequent analyses for forward trimmers were made with an assumed trimmer 
hinge margin of 20 percent, which would place the hinge axis at about 5 percent 
of the trimmer mean aerodynamic chord. 

Table IX lists the results for two modified versions of Configuration 
1A, with moment arms of 1.5 and 1.0 chord lengths. The free-trimmer hinge 
margin is 20 percent in both cases. 
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TABLE IX. RESPONSES OF MODIFIED FORWARD TRIMMER CONFIGURATIONS 


Trimmer 


Response 

ratio 





moment 

arm, 

Rms 

compared 
to fixed 

Mode A 


Mode B 


wing 

chords 

turbulence 

response 

wing, 

percent 

Roots 

Damping 

ratio 

Roots 

Damping 

ratio 

1.5 

0.0178 

98 

-0.431 + j 3.35 

0.131 

-3.09 + j 7.03 

0.402 

1.0 

0.0154 

81 

-0.445 + j 3.29 

0.134 

-5.63 + j 10.67 

0.467 


As a basis of comparison, the rms gust response for the pure free 
wing of Configuration IB is 0.00989, or just 52 percent of the fixed-wing 
response. Even the smaller moment arm case falls far short of the basic free 
wing, from a gust alleviation standpoint, although a reduction of 19 percent 
from the fixed-wing response is significant. 

It should be noted that the two oscillatory modes become more widely 
separated as the pitching inertia of the wing assembly decreases. With de- 
creasing moment arm, the frequency of mode B increases, while mode A changes 
only slightly. 

If the fuselage assembly is considered by itself, free to heave and 
pitch, its characteristic motion is defined by a cubic equation with a well- 
damped real root and an oscillatory mode with roots of -0.712 + j 3.43 at this 
flight condition. It appears that mode A of the complete aircraft is tending 
toward this isolated fuselage oscillatory mode; but despite the numerical 
separation of modes A and B, both modes are evident in the load factor response 

history shown in figure 6 for the one chord length moment arm case. 

The handling qualities specification of MIL-F-8785B (ASG) is based 
upon the assumption that the short-period response is essentially that of a 
simple second order system. In the case at hand, the short-period motion is 
clearly the superposition of two separate oscillatory modes, raising a question 

of the direct application of the simplified standards to the f ree-wing/free- 

trimmer aircraft. 

Aft-trimmer configurations. - The most powerful aft-trimmer config- 
uration examined was one based upon Conf iguration 1A, but with the trimming 
surface mounted two wing chord lengths aft of the wing hinge axis. 
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FIGURE 6. RESPONSE TO STEP CONTROL INPUT OF MODIFIED ■ 

CONFIGURATION 1A AIRCRAFT WITH FORWARD , 

TRIMMER MOMENT ARM OF ONE WING CHORD LENGTH. I 

i 

The adverse short-period damping of this configuration has been t 

mentioned previously under the discussion of longitudinal characteristic modes, 
and a time history of the response to a step control input is shown in 
figure 4. 

The damping ratio of the dominant mode in figure 4 is only 0.027, 
which is obviously well below the standards of MIL-F-8785B (ASG) , which requires 
at least 0.15 for level 3 flying qualities. The rms responses to continuous 
turbulence were not obtained for this configuration because of excessive com- 
puting time caused by the large spike in the power spectrum of the normal load 
factor response. 

To improve the behavior of the aft-trimmer configurations, smaller 
moment arms were explored. Unlike the forward trimmer cases, however, the 
trimmer hinge margin was maintained at the 6 percent value to retain a rela- 
tively high value of maximum trimmer lift coefficient. Furthermore, the 
trimmers in all subsequent cases were assumed to be mounted on the wingtips, 
as seen in Configuration 1C (fig. 2(c)). 

Aside from esthetic improvement, the tip-mounted trimmers are 
superior, from a gust alleviation standpoint, to surfaces mounted behind the 
wing. The aft trimmer experiences a transient lift change, upon gust penetra- 
tion, just like the forward surface of Configuration 1A; but in this case the 
transient force on the aft trimmer is beneficial, since it produces a favorable 
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pitching moment on the main wing. With the outboard location, the trimmer is 
exposed to a strong upwash increment from increased circulation in the trailing 
vortex system of the wing. The upwash increases the magnitude of the transient 
trimmer force and enhances the relieving pitching moment on the wing. If the 
trimmer is mounted directly behind the wing, the influence of the wing appears 
as an increase in downwash (for an upward gust), and the alleviating tran- 
sient moment from the trimmer is reduced. 

Table X lists the behavior of aircraft in cruise flight with tip- 
mounted aft trimmers for two moment arm lengths and for an area ratio of 1/6, 
as before. 



TABLE X. 

BEHAVIOR OF 

TIP-MOUNTED AFT- 

-TRIMMER 

CONFIGURATIONS 


Moment 

arm, 

Rms 

Response 
compared 
to fixed 

Mode A 


Mode B 


wing 

chords 

turbulence 

response 

wing, 

percent 

Roots 

Damping 

ratio 

Roots 

Damping 

ratio 

-1.5 

0.0132 

69 

-0.509 + j 4.54 

0.111 

-1.77 + j 2.97 

0.512 

-1.0 

0.0104 

54 

-1.11 + j 3.32 

0.317 

-1.67 + j 5.28 

0.302 


The gust alleviation of the one chord length moment arm configu- 
ration (Configuration 1C, fig. 2(c)) is excellent, inasmuch as it is virtually 
identical to the pure free wing, which reduced the response to 52 percent of 
the fixed-wing value. Furthermore, the computed damping ratios of both modes 
exceed the level 3 requirements (0.15) of MIL-F-8785B (ASG) , and almost meet 
the level 1 standard (0.35). 

The larger moment arm length (one and a half wing chords) is infe- 
rior to the smaller moment arm case with respect to both turbulence response 
and short-period damping. 

Since the literal application of the MIL-F-8785B (ASG) short-period 
damping standards is open to question because of the bimodal nature of the re- 
sponse, time histories of load factor response to step control displacement 
were computed and are shown in figure 7. The response of the smaller moment 
arm configuration, in particular, would appear to be acceptable from an intui- 
tive viewpoint . 
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FIGURE 7. COMPARISON OF LOAD FACTOR RESPONSES TO STEP CONTROL INPUT 
FOR TWO CONFIGURATIONS WITH TIP -MOUNTED AFT TRIMMERS 

Effect of configuration on maximum lift coefficient. - The principal 
advantage of the external trimmer concept over the basic free wing is the fea- 
sibility of improving the maximum trimmed lift coefficient of the aircraft. If 
this were the only consideration, the forward trimmer configuration with a 
large moment arm would be attractive, because it would permit the use of rela- 
tively powerful wing flaps and would contribute an upward force to augment the 
wing lift. 

Unfortunately, very powerful external trimming surfaces, mounted 
either forward or aft of the wing, have been found to have adverse dynamic 
characteristics within the scope of this study, as previously discussed. 

These circumstances prompted a brief investigation of the methods of increasing 
the trimming power with relatively modest moment arms. 

The fundamental problem is the large negative pitching moment, which 
inevitably accompanies the trailing-edge flap deflection required for high wing 
lift coefficients. The wing pitching moment about its hinge axis is given by 
the following equation: 

Cm binge ' C %.c. ' <™> C K, <2) 
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C m is always large and negative for large trailing-edge flap 

a • c • 

deflections; for example, it has a value of about -0.40 for 20 percent flaps 
deflected 40°. 

The moment available from the external trimming surface to counter 

the wing moment of equation (2) is as follows: 

X’ S 

C = _c / 

m trimmer - S 1 
c c 

Since the maximum lift coefficient of a free trimming surface con- 
trolled by a trailing-edge tab can be expected to be about 1.0, the maximum 
trimmer pitching moments which can be obtained are given in table XI. 


TABLE XI. MAXIMUM TRIM MOMENT CAPABILITY 
FOR PLAIN TRIMMERS 


K 

c 


!c 

S 

Cm . 
“trimmer 

1.0 

1/6 

0.167 

1.5 

1/6 

0.250 

2.0 

1/6 

0.333 

1.0 

1/5 

0.200 

1.5 

1/5 

0.300 

2.0 

1/5 

0.400 


It should be noted that only the last case listed in table XI could 
provide sufficient power to trim the constant pitching moment (C m ) of a 
wing with 20 percent flaps deflected 40°, should large deflection be required. 
In that case, no excess moment would be available to counter the moment due to 
lift in equation (2) . 

In reference 2, the lift capability of a pure free-wing aircraft was 
improved by using a leading-edge Handley-Page slat on the NACA 23012 airfoil 
section, using data from reference 3. This device yielded a computed maximum 
trimmed of 1.46 for a rectangular wing of aspect ratio 6. If such a device 
were used on the free-trimmer surfaces, the increased trimmer power would pro- 
duce the maximum moments given in table XII. 

Although the trimming power is increased significantly by the trimmer 
leading-edge slat (46 percent) , large wing flap deflections still cannot be 
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TABLE XII. MAXIMUM TRIM MOMENT CAPABILITY FOR TRIMMERS 
EQUIPPED WITH HANDLEY-PAGE LEADING-EDGE SLAT 


f 

X c 

c 

S c 

S 

Q 

^trimmer 

1.0 

1/6 

0.243 

1.5 

1/6 

0.365 

2.0 

1/5 

0.487 

1.0 

1/5 

0.292 

1.5 

1/5 

0.438 

2.0 

1/5 

0.584 


tolerated for moment arms of 1.5 chord lengths or less. For this reason, only 
smaller flap deflections are considered in the following calculations. 

To provide parametric information, the data in figure 119 of refer- 
ence 4 were used to derive a simple, approximate relationship between the 
maximum section lift and pitching moment coefficients of the NACA 23012 air- 
foil with variation in the deflection of a 20 percent slotted flap. This 
relationship is: 


' J, m j n 


= 0.35 - 0.24 Cl 


(C 1 >1.5) 


(4) 


max 


u a . c . J -max 

If it is assumed that the finite wing lift coefficient is 80 percent 
of the maximum section lift coefficient, equation (4) becomes 


C m = 0.35 - 0.30 C T 
m a . c . L max w 


( 5 ) 


If the flap is deflected just sufficiently to produce Cl 


max. 


the 


substitution of equation (5) into equation (2) provides a convenient expression 

for the total pitching moment of the wing about its hinge axis at Cl • For 

maxr 

the nominal hinge margin of 0.06, the expression is 


w 


C m , . = 0.35 - 0.36 C T 0 

m hinge L max w 

Equating the hinge moment of equation (6) to the available trimmer 
moments of tables XI and XII provides an estimate of the maximum permissible 
wing lift coefficient. Ignoring the contribution of the fuselage assembly, 
which is a function of fuselage center of gravity and trimmed attitude, the 
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total trimmed airplane lift coefficient is 

s c 

Ct — Ct 4- —pr~ Ct ( 

L total Her — S L c 

where the plus sign is used for the forward trimmer configurations and the 
minus sign is used for aft trimmers. 

Table XIII lists the total trimmed airplane lift coefficients ob- 
tainable with the 20 percent wing flaps deflected just enough to provide the 
usable maximum wing lift coefficient. 

TABLE XIII. MAXIMUM AIRPLANE LIFT CAPABILITIES FOR WING 
EQUIPPED WITH 20 PERCENT SLOTTED FLAPS 


*'c 

s 

c 

Usable 

Cr 

lj max w 

. 5 f 

degrees 

Total 

c L 

c 

s 

Forward 

trimmer 


Aft trimmer 



[Plain trimmer] 




1.0 

1/6 

1.44 

6.6 

1.60 


1.27 

1.5 

1/6 

1.67 

13.0 

1.83 


1.50 

2.0 

1/6 

1.90 

19.4 

2.06 


1.73 

1.0 

1/5 

1.53 

9.1 

1.73 


1.33 

1.5 

1/5 

1.81 

16.8 

2.01 


1.61 

2.0 

1/5 

2.08 

24.5 

2.28 


1.88 



[Leading-edge slat on trimmer] 



1.0 

1/6 

1.65 

12.4 

1.89 


1.41 

1.5 

1/6 

1.99 

21.8 

2.23 


1.75 

2.0 

1/6 

2.33 

31.3 

2.57 


2.09 

1.0 

1/5 

1.78 

16.2 

2.07 


1.49 

1.5 

1/5 

2.19 

27.5 

2.48 


1.90 

2.0 

1/5 

2.59 

38.7 

2.88 


2.30 


Another possible approach to a high-lift capability would 

be to 

provide 

a leading-edge slat on the 

wing itself, in combination 

with 

trailing- 

edge flaps. 







Reference 

3 provides data for the NACA 

23012 airfoil 

with 

a Handley- 

Page leading-edge slat and a 26 percent slotted flap. From the data in 

figure " 

1 of reference 3, it may be 

deduced that the aerodynamii 

c center of the 
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airfoil section shifts forward to about the 19 percent chord position as the 1 
leading-edge slat is extended. Since this is also the position of the nominal 
hinge axis location, it would be necessary to move the hinge forward to provide 

t 

static wing stability with the slat extended. If a new hinge location is 
selected at 15 percent chord, the hinge margin would be 10 percent with the 
slat retracted and 4 percent with the slat extended. I 

Using the data in reference 3, and again assuming that the maximum 
wing lift coefficient is 80 percent of the maximum section value, the relation- 
ship between the pure pitching moment about the aerodynamic center (0.19c) and 
the maximum lift coefficient (both of which are functions of flap deflection) 
is, with the leading-edge slat extended, 

C = 0.665 - 0.438 C T (8) ! 

m L i 

a . c . max 

w j 

If the new hinge margin is 0.04, * 

C = C - 0.04 C (9) 1 

m, . m L 

range a.c. w : 

Substituting equation (8) into equation (9) , the hinge moment about ! 
the wing hinge when the wing is at the maximum lift coefficient as determined 
by flap deflection, is 

C = 0.665 - 0.478 C T (10) 

m, . L 

hinge max 

w 

Equating the hinge moment of equation (10) with the available 
trimmer power in tables XI and XII yields the data listed in table XIV. 

The data for the various high-lift configurations are summarized in 
figure 8 for the smaller of the two area ratios considered. It may be seen 
that the maximum lift coefficient of the pure free wing is easily exceeded, 
even for aft trimmers with relatively small moment arms. 


i 


Weight penalties. - In an attempt to assess the weight penalties for 
each of the free-wing /f ree-trimmer combinations studied, each wing-trimmer com- 
bination was added to one basic fuselage incorporating all fixed and disposable 
weight items for both the two seat light aircraft and the single seat agricul- 
tural aircraft (tables III and VI) . 

For the two seat light aircraft, given a fixed-wing baseline (Con- 
figuration ID in table III), the following rules were used: all non-wing/ 
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TABLE XIV. MAXIMUM AIRPLANE LIFT CAPABILITIES FOR WING EQUIPPED 
WITH LEADING-EDGE SLAT AND FLAPS 


Usable Total 


A 

c 

c 

Sc 

S 

L 

max 

w 

6F 

degrees 

Forward 

trimmer 

Aft trimmer 




[Plain trimmer ] 



1.0 

1/6 

1.90 

6.9 

2.07 

1.74 

1.5 

1/6 

2.09 

12.8 

2.26 

1.92 

2.0 

1/6 

2.28 

18.8 

2.45 

2.11 

1.0 

1/5 

1.98 

9.3 

2.18 

1.78 

1.5 

1/5 

2.21 

16.4 

2.41 

2.01 

2.0 

1/5 

2.43 

23.6 

2.63 

2.23 



[Leading-edge slat on trimmer] 


1.0 

1/6 

2.08 

12.3 

2.32 

1.83 

1.5 

1/6 

2.35 

21.1 

2.60 

2.11 

1.0 

1/5 

2.19 

15.8 

2.48 

1.90 

1.5 

1/5 

2.52 

26.3 

2.81 

2.23 

trimmer 

weight items 

were to remain constant , wing 

area was to 

remain constant 


at 10.1 square meters (108.6 square feet), trimmer area was to remain constant 
at 1.7 square meters (18.1 square feet), wing/boom and trimmer free surfaces 
were to be externally counterbalanced, and useful load was to be constant. The 
free-wing/forward-trimmer configuration (configuration 1A in fig, 2(a)) had 
an empty weight 26 percent greater than its fixed-wing counterpart and 18 per- 
cent greater than a comparable pure free-wing aircraft. The f ree-wing/af t 
trimmer conf iguration (configuration 1C in fig. 2(c)) had an empty weight 
32 percent greater than its fixed-wing counterpart and 23 percent greater than 
a comparable pure free wing. (These weights reflect no attempt to minimize 
the weight gain of various configurations by making efficient use of some items 
of fixed weight as counterweights for wings and trimmers.) These weights 
essentially reflect the penalties incurred in modifying a fixed- or pure free- 
wing aircraft to free-wing/ free-trimmer configurations without redesigning the 
basic aircraft. 
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Moment Arm from Wing Hinge, wing chords 

(a) Forward trimmers. (b) Aft trimmers. 

FIGURE 8. MAXIMUM TRIMMED AIRPLANE LIFT COEFFICIENTS OF VARIOUS 
FREE- WING/ FREE-TRIMMER CONFIGURATIONS . 


Following the same rules for the single seat agricultural aircraft, 
the results are somewhat different. Considering the fixed-wing configuration 
(configuration 2D in table IV) as a baseline, the free-wihg/ aft- trimmer con- 
figuration (configuration 2C in fig. 3(c)) had an empty weight 21 percent 
greater than configuration 1A, its fixed-wing counterpart, and 7 percent 
greater than the comparable pure free-wing aircraft. The f ree-wing/f orward- 
trimmer configuration (configuration 2A of fig. 3(a)) had an empty weight 
8 percent greater than its fixed-wing counterpart but 5 percent less than a 
comparable pure free-wing aircraft. 

This comparison brings up important design considerations. First, 
an aft-trimmer configuration incurs the highest weight penalty because the 
center of gravity of the trimmer surface is aft of the wing hinge line and 
does nothing to alleviate the counterweight needed to balance the free wing 
about its hinge line. Secondly, a forward- trimmer configuration is more 
efficient from a weight standpoint in that it incurs a lower weight penalty 
than an aft trimmer. Since the weight of the forward trimmer offsets the 
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weight needed to balance the wing about its hinge point, it can, if properly 
sized and placed, provide a lighter total wing weight than a pure free wing. 

Obviously, a tradeoff must be made on a conf iguration-by-conf ig- 
uration basis to size the external trimming surface, and to decide whether 
one is, in fact, necessary. For static balance, the boom/trimmer combination 
must create just enough moment about the wing hinge axis to bring the wing/ 
boom/ trimmer center of gravity to the wing hinge axis. On the other hand, the 
wing/boom/trimmer pitch inertia must be as low as possible to provide accep- 
table handling qualities. The trimmer must also be close enough to the wing to 
provide a ride quality comparable to that of pure free-wing aircraft and far 
enough ahead of the wing to provide adequate control power for generating 
high trimmed wing lift coefficients. 


Sensitivity to Mass Imbalance 

All of the configurations examined to this point have incorporated 
the longitudinal mass balancing of the three major components. Although the 
fuselage center of gravity is calculated to be approximately 0.61 meter (2 feet) 
below the wing hinge, the longitudinal displacement with respect to the hinge 
axis is zero. The centers of gravity of the wing/boom/trimmer assembly and 
the trimming surface coincide exactly with their respective hinge axes. 

The effects of mass imbalance in each of the three major components 
were examined by computing the loci of the characteristic roots as the centers 
of mass deviated in the fore and aft directions. For consistency, the aircraft 
in Configuration 1C was used throughout, with a tip-mounted aft trimmer posi- 
tioned one wing chord length behind the wing hinge axis. 

Two flight conditions were examined for each case: the nominal 
cruise condition, which was flight at 241.4 kilometers per hour (150 miles 
per hour) at an altitude of 1829 meters (6000 feet), and an approach condition, 
which was flight at 121.0 kilometers per hour (75 miles per hour) at sea level. 
The assumed approach speed is 30 percent above the minimum speed obtainable 
with a slatted trimmer and wing flaps and slats (case 4 in fig. 8(b)) at a 
gross weight of 810.88 kilograms (1788 pounds). For consistency with the high- 
lift configuration, the wing hinge margin is assumed to be 4 percent during 
approach, as discussed in the preceding section. The combination of low 
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airspeed and reduced wing hinge margin is intended to provide the highest 
sensitivity to mass imbalance because of the higher ratio of inertial forces 
to aerodynamic forces. 

Fuselage center of gravity location. - While the wing assembly and 
trimming surface can be designed to have an invariant center of gravity, the 
fuselage assembly contains all variable and disposable load items, and provi- 
sion must be made for a reasonable center of gravity range. 

In cruise, the aircraft of Configuration 1C could tolerate a 
forward displacement of about 0.6 meter (2 feet), which is 47 percent of the 
mean aerodynamic chord, and an aft movement of even greater magnitude, accord- 
ing to the root locus plots in figure 9. At the forward limit, mode A damping 
deteriorates to the level 3 standards of MIL-F-8785B (ASG) , and the phugoid 
mode becomes nearly neutrally stable. As the fuselage center of gravity moves 
aft from the wing hinge, the only adverse effect appears to be some degrada- 
tion of the damping of mode B. 

In the approach case, the effects are more pronounced, and the 
movement of the characteristic roots, as seen in figure 10, has little resem- 
blance to the cruise condition. The damping ratio of mode A is slightly 
inferior to level 3 standards even for a mass balanced fuselage, and even a 
small aft movement is sufficient to cause a totally unacceptable dynamic 
instability. In addition, a forward displacement of less than 0.3 meters 
(1 foot) would be sufficient to produce an unstable phugoid motion. 

As mentioned earlier, the coupling of fuselage and wing motion is 
much more in evidence with the external trimmer arrangements than with the pure 
free wing concept. For this reason, and because Configuration 1C appeared to 
be only marginally acceptable in approach, a brief excursion was made to 
explore the effects of changes in the aerodynamics of the fuselage assembly. 
Specifically, a large-tail version of Configuration 1C was hypothesized in 
which the fuselage stabilizer moment arm was held constant but the effective 
area of the stabilizer was doubled. This modification substantially increased 
the static angle of attack stability of the fuselage assembly, the fuselage 
pitch damping coefficient, and the slope of the fuselage lift curve. 

The variation of fuselage center of gravity was then repeated for 
the large-tail version. The resulting root loci are plotted in figure 11 for 
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MIL - F- 8785B ( ASG ) 
minimum damping ratio 
requirements 












MIL - f - 8785B ( ASG) 



FIGURE 11. EFFECT OF FUSELAGE CENTER OF GRAVITY ON LARGE- 
TAIL VERSION OF CONFIGURATION 1C. APPROACH 
CONDITION WITH 4 PERCENT HINGE MARGIN. 
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approach and figure 12 for cruise. Again, the approach case is the more sensi- 
tive, but the fuselage center of gravity could deviate 0.3 meter (1 foot) for- 
ward and 0.6 meter (2 feet) aft before violating the level 1 phugoid standards 
or the level 3 short-period damping standards, respectively. These deviations 
would correspond to a fuselage center of gravity range of about 71 percent of 
the wing mean aerodynamic chord. 

Although the approach dynamics of the large-tail version are im- 
proved and the sensitivity to fuselage imbalance is decreased, the modification 
would not be without its price. Comparing figures 9 and 12, the damping ratio 
of mode B in cruise flight is adversely affected; and more importantly, the 
turbulence response is inferior to the original configuration, Configuration 1C 
The turbulence response will be covered in a subsequent section. 

Because of the coupling between the dynamics of the fuselage and 
wing assembly with external trimmer, it is difficult to come to definite con- 
clusions regarding the effects of fuselage center of gravity displacement. 
Nevertheless, from the four data points provided by the two flight conditions 
and two tail sizes, it may be concluded that forward displacement of the fuse- 
lage center of gravity decreases the damping of the phugoid mode and that aft 
displacement decreases the damping of one of the short-period modes, either 
A or B. Further, the effect of fuselage imbalance is more pronounced at low 
speeds and depends upon the aerodynamic design of the fuselage assembly. Sig- 
nificantly, however, the permissible range can be comparable to fixed-wing 
aircraft. 

Wing assembly imbalance. - The center of gravity of the complete 
wing assembly, including the trimmer mass, was permitted to move through a 
distance of 15 percent of the wing chord fore and aft of the wing hinge. The 
resulting root loci are plotted in figure 13 for the cruise condition and 
figure 14 for the approach condition. 

As shown by figure 13, the prescribed center of gravity displace- 
ment has no significant effect on any of the characteristic modes. For the 
approach condition, on the other hand, large forward center of gravity displace 
ments (in excess of about 12 percent of wing chord) cause an instability in 
mode A. As plotted in figure 14, no serious detrimental effects are noted for 
aft displacements . 
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FIGURE 12. EFFECT OF FUSELAGE CENTER OF GRAVITY ON LARGE-TAIL 
VERSION OF CONFIGURATION 1C. CRUISE CONDITION WITH 
10 PERCENT HINGE MARGIN. 
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FIGURE 13. EFFECT OF WING /BOOM/ TRIMMER IMBALANCE ON CHARACTER- 
ISTIC ROOTS OF CONFIGURATION 1C. CRUISE FLIGHT WITH 
6 PERCENT HINGE MARGIN. 










(b) Phugoid mode. 


FIGURE 14. EFFECT OF WING/ BOOM/ TRIMMER IMBALANCE ON CHARACTERISTIC 
ROOTS OF CONFIGURATION 1C. APPROACH CONDITION WITH 4 
PERCENT HINGE MARGIN. 
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Although the configuration under examination appears to be very 
tolerant to wing center of gravity displacement, the adverse effect of aft 
movement of the center of gravity during landing was not analyzed. In ref- 
erence 2, it was concluded that an aft wing panel center of gravity would make 
a smooth landing very difficult because of the tendency for the free wing to 
pitch upward in response to a finite landing impact. The abrupt increase in 
angle of attack, in turn, would cause a bounce which could lead to a series 
of encounters with the runway before sufficient airspeed loss was achieved to 
remain on the ground. The landing impact simulation was not repeated for the 
f ree-wing/free-trimmer concept, but there is no reason to expect the results 
to be substantially different. 

It may be concluded that small variations in wing center of gravity 
(of the order of a few percent of wing chord) will have no significant effect 
on the in-flight characteristic modes, but aft center of gravity locations 
should be avoided to facilitate smooth landings. 

Trimmer imbalance. - The trimmer center of gravity was displaced 
up to 5 percent of the trimmer chord fore and aft of its hinge axis, and the 
characteristic modes were found to be very sensitive to this parameter. 

An aft imbalance of the trimmer is especially to be avoided, as 
evidenced by the movement of the phugoid roots in figure 15 for cruise, and in 
the instability of mode A in figure 16 for the approach condition. For the 
latter case, an aft displacement of only 0.4 percent of the trimmer chord pro- 
duces a dynamic instability. For forward center of gravity movement, greater 
tolerance is exhibited before the phugoid mode is destabilized for the ap- 
proach condition. 

Since the damping ratio of mode A in the approach condition is only 
marginally acceptable, even with no imbalance, further numerical experimenta- 
tion was performed to insure that this significant result was not excessively 
configuration dependent. Again, the large-tail version of Configuration 1C was 
used. In this case, mode A improved but mode B became unstable, with an aft 
displacement only slightly in excess of 1 percent of the trimmer chord. 

Because of the criticality of this parameter, the numerical experi- 
ment was repeated with a forward trimmer configuration to explore possible 
differences. Configuration 1A was selected with a forward surface at a moment 
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arm of two wing chords. With regard to the consequences of aft trimmer im- 
balance, the results were virtually the same, a dynamic instability of mode B 
occurring with an aft imbalance of less than 1 percent of the trimmer chord. 

It would appear that the effect of aft imbalance is the same 
whether the trimmer is located fore or aft of the wing hinge. A reasonable 
explanation for this observation is that, in both configurations, an aft im- 
balance of the panel tends to deflect the trimmer in the direction that en- 
hances the angular acceleration of the wing. 

In conclusion, no aft imbalance of the trimmer surface should be 
permitted. A small forward imbalance (of the order of one percent of trimmer 
chord) may be tolerated, although the effect may be felt in reduced damping of 
the phugoid at low speeds. 

General Comments on Configuration Options 

Although the forward trimmer configurations afford higher lift 
capabilities (fig. 8), the aft-trimmer configurations appear to offer inher- 
ently better gust alleviation. Furthermore, for the stick-fixed longitudinal 
motions analyzed in this study, the aft-trimmer configurations exhibit better 
dynamic behavior (compare fig. 6 and 7). 

Figure 17 summarizes the normal load factor responses to continuous 
turbulence for both types of configurations for the cruise condition. Config- 
uration 1C, with a tip-mounted aft trimmer with a moment arm of one chord 
length, was selected for the analysis of imbalance effects because it appears 
to be representative of promising configuration choices. For completeness, 
the turbulence responses were also calculated for the approach condition for 
this configuration; these responses are shown in figure 18. 

The adverse effect of the free-wing principle on fuselage pitching 
response is consistent with the observations of references 1 and 2. The devel- 
opment of passive or active fuselage pitch dampers would be an interesting 
technical challenge because of the influence of fuselage characteristics on 
total aircraft behavior which has been noted for the free-wing/ free- trimmer 
conf igurations . 

One disadvantage of aft trimming surfaces is the necessity for 
greater ballast forward of the wing hinge to achieve mass balance of the wing 
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FIGURE 17. COMPARISON OF RESPONSES TO CONTINUOUS TURBULENCE. 



FIGURE 18. 


COMPARISON OF RESPONSES OF FIXED-WING 
AIRCRAFT AND TWO VERSIONS OF 
CONFIGURATION 1C. 


assembly. It seems possible, however, to reduce the severity of this penalty 
by clever design of specific aircraft. For example, landing lights, avionics 
equipment, batteries, and other constant-mass fixed items of equipment could, 
in principle, be placed in the forward wingtip booms to reduce the required 
mass of additional inert weight. In addition, judicious choices of wing sweep 
and taper may move the wing center of gravity forward with respect to the wing 
hinge axis. 

Another intriguing possibility for tip-mounted aft trimmers is the 
use of these surfaces for lateral control. Although no lateral-directional 
analysis was included in this study, the trimming surfaces, moved differen- 
tially, could be expected to provide powerful roll control. 


CONCLUSIONS 


The following conclusions may be drawn from this investigation: 

(1) For the trimmer area ratio considered (1/6), the most 
promising configuration employs wingtip-mounted trimming 
surfaces placed aft of the wing hinge line with a moment 

arm of one wing chord length. Of the configurations examined 
in this study, this arrangement alone could provide excell- 
ent alleviation of vertical gust loads while exceeding the 
maximum lift capability of pure free-wing configurations, 
and while meeting fundamental criteria for the stability 
of the stick-fixed longitudinal modes. 

(2) For vertical gust alleviation, forward trimmers are inferior 
to aft-mounted surfaces because of adverse wing pitching 
moments caused by transient aerodynamic forces on the trim- 
ming surfaces. 

(3) Mass balancing of the trimmer surface about its hinge axis 
is vital for precluding adverse effects on the stability of 
the characteristic modes. In particular, aft imbalance 
must be avoided. 

(4) Longitudinal displacement of the center of gravity of the 
fuselage assembly appears to be more significant for free- 
wing/f ree-trimmer configurations than for pure free-wing 
aircraft. Forward displacement decreases the damping of the 
phugoid mode while aft displacement decreases the damping 

of one of the short-period modes. The effect of fuselage 
imbalance is more pronounced for slow-speed flight, and the 
sensitivity depends upon the aerodynamic design of the fuse- 
lage assembly. 
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(5) Small variations in the wing assembly center of gravity (of 
the order of a few percent of wing chord) have no signifi- 
cant effect on the in-flight characteristic modes, but center 
of gravity locations aft of the wing hinge axis should be 
avoided to facilitate smooth landings. 

(6) Forward-trimmer configurations are more efficient from a weight 
standpoint than aft trimmers, and could, if properly sized and 
placed, provide a lighter total wing weight than a pure free- 
wing. The aft-trimmer configuration incurs a higher weight 
penalty because of the additional counterweight needed to 
balance the wing assembly about its hinge axis. 


APPENDIX A. DEVELOPMENT OF EQUATIONS OF MOTION 


Introduction 

In deriving the equations of motion, each free surface and the fuse- 
lage assembly are considered free bodies at first. After the individual sets of 
equations are written with respect to the most convenient axis systems, they are 
combined into a single set, referred to the wing hinge axis system. The 
equations are consolidated by eliminating the common forces and moments that act 
between the various components. The equations are then linearized for conven- 
ience in the analysis. 


Symbols 


Symbols that are defined each time they are used have been omitted 
from this list. 


a 

F 

F 

g 

H 

m 

M 


2 2 

linear acceleration, meters/second (feet/second ) 
force, newtons (pounds) 

reversed effective force, newtons (pounds) 

2 2 

acceleration due to gravity, meters/second (feet/second ) 
angular momentum, radians/second 

inertia tensor (shown in equation (31) for forward trimmer and 
after equation (44) for wing-boom assembly), kilogram-meters^ 
(slug-feet^) 

moment , newton-meters (foot-pounds) 


M 

reversed effective moment, newton-meters (foot-pounds) 

m 

mass, kilograms (slugs) 


roll rates about x,y,z axes, respectively, radians/second 

R 

2 2 

inertial linear acceleration, meters/second (feet/second 

u,v,w 

velocity vector components about x,y,z axes, respectively, 
meters /second (feet/ second) 

V 

velocity, meters/second (feet/second) 

X, Y,Z 

stability axes 

X H’ y H’ Z H 

wing hinge axes 

x' ,y' ,z' 

wing panel axes 

yi y 1 7 i 

CH’ CH’ CH 

trimmer hinge axes 
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X",Y",Z" 

6 

P 


1 ,1 ,1 

x y z 

Subscripts and 


c 

eg 


CH 

F 

H 

0 

O f 


P 

W 

WB 

WH 

x,y,z 

i 

?t 


trimmer panel axes 

deflection angle of surface, radians 
position vector, meters (feet) 
angular velocity, radians/second 
unit vectors along x,y,z axes, respectively. 

superscripts : 

free trimmer surface 
center of gravity 
trimmer hinge 
fuselage 
wing hinge 

wing hinge line intersection with aircraft plane of symmetry 

trimmer hinge line intersection with aircraft plane of 
symmetry 

wing panel 

wing 

wing/boom assembly 
wing hinge 

axes other than stability 
wing panel axis system 
trimmer axis. 


Coordinate Systems 

The six coordinate systems described below were employed. 

Conventional stability axis system . Following standard practice, the 
basic set of coordinates for describing the aircraft motion has its origin at 
the center of gravity of the complete aircraft. The x-axis is aligned with the 
velocity vector of the aircraft in the reference condition, the y-axis extends 
to the right of the plane of symmetry, and the z-axis completes the right-hand 
set. These coordinates are fixed in the aircraft and rotate with it. The 
orientation of the stability axis system with respect to an inertially fixed 
reference is defined by three standard Euler angles. The sequence of rotation 
used to define these angles is: (1) rotation about the z-axis through the yaw 

angle, (2) rotation about the y-axis through the pitch angle, 0, and (3) 

rotation about the x-axis through the roll angle, <f>. 
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Hinge axis system . The hinge system of axes, haS its ori S in 

in the plane of symmetry of the aircraft. The positive y -axis coincides with 
the axis of rotation of the right wing panel. For simplicity, the wing panels 
are assumed to have no geometric dihedral. Consequently, z lies in the plane 
of symmetry, and the negative y^-axis coincides with the axis of rotation of the 
left wing panel. The hinge axis system is parallel to the stability axis system., 
and is therefore fixed in the fuselage assembly for a given flight condition. 

Panel axis system . The panel axis system x 1 , y f , z ' , is similar to 
the hinge axis system, but rotates with the wing panels. When dealing with the 
wing panels, the panel axis system is rotated about the y^-axis through the dis- 
placement angle 6 p . 

Forward trimmer hinge axis system . The forward trimmer hinge axis 

system, X* Y' Z' is fixed with respect to the wing panel axis system and 
CH CH CH 

varies from the stability axis system by 6 p only. Its origin is at the inter- 
section of the forward trimmer hinge line and the aircraft plane of symmetry. 

Forward trimmer panel axis system . The forward trimmer system of 
axes X M , Y" > Z M has its origin at the same point as the trimmer hinge axis 
system. The forward trimmer panel axis system, however, is fixed in the forward 
trimmer panel, its x-axis being the forward trimmer chord line at the aircraft 
plane of symmetry. This axis system is free to rotate from the forward trimmer 
hinge system through the displacement angle 6^ , only. The Y n - and Y^-axes 
coincide and both are parallel to the y'-, y^, y p -, and Y-axes at all times. 

Fuselage axis system . The fuselage axis system x p ,y p ,z p is fixed 

within the fuselage at the wing hinge. Its origin is at the intersection of 

the aircraft plane of symmetry and the wing hinge line, since this is the point 

about which the fuselage rotates. The x-axis is parallel to the fuselage 

r 

reference line (FRL) , positive forward. The axis y p is coincident with y p and 
y T , and x^ and z vary from x and z by the fuselage pitch angle, e. The fuse- 
lage pitch angle is usually small enough to be neglected. To assess the effect 
of fuselage pitch angle on aircraft characteristics, e was initially set at a 
value other than zero and fuselage aerodynamic and inertial characteristics 
recalculated to reflect its static value. 
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Whenever equations of motion are to be formulated, they always derive 
from the application of two principles; namely, the principle of motion of the 
mass center and the principle of moment of momentum. The principle of motion of 
the mass center applies only to centers of gravity and this was carefully 
applied to each of the free bodies (fuselage, wing/ boom and trimmer). The 
principle of moment of momentum is valid for points of rotation other than the 
mass center, if compensation is made for the offset center of gravity. This 
compensation was generally necessary for each free body, since the various 
auxiliary axis systems were not generally located at centers of gravity. 

Unlike the case for a conventional aircraft, with the free-wing air- 
craft the equations of motion cannot be summarily referenced to the stability 
axes centered at the aircraft center of gravity. In fact, the location of the 
all-up aircraft center of gravity is variable, dependent upon the relative 
position of the wing/boom/ trimmer . This variability has not been ignored in the 
current analysis. The hinge axis system can be used as the "principal 
reference" rather than the stability axis system. It is fixed with respect to 
the fuselage and parallel to the stability axis. All force and moment contri- 
butions from the various free bodies can be referenced/coordinated with respect 
to this axis system, even though it is not located at the aircraft center of 
gravity, nor at the center of gravity of any of the various free bodies. 

Free-Body Equations 

General procedure . To derive forward trimmer and wing panel force 
equations, first the position vector, p, from the panel hinge axis (coordinate 
system origin) to the panel center of gravity is written in the most convenient 
coordinate system. This allows the vector to be fixed. The position vector is 
then transferred to the wing hinge axis system and differentiated twice with 
respect to time. To this is then added the time differential of inertial 
velocity. 

To derive forward trimmer and wing moment equations, first the posi- 
tion vector from the panel hinge axis to the panel center of gravity is written 
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I 


in the appropriate axis system, as before. Then the momentum term, mR, and the 
angular momentum term, H, are written in the appropriate axis systems. The 
I vector p, and the terms mR and H are transferred to the wing hinge axis system 
before being differentiated with respect to time as necessary to write the 
I appropriate moment equations. 

The above steps also apply to the fuselage moment equations, but these 
equations involve only the fuselage and stability axis systems. 

The aircraft will be broken down into three free bodies for analysis: 

| 

the forward trimmer, wing/boom assembly and fuselage. 

Forward trimmer panel force equations . In the forward trimmer panel 

| axis system in figure 19, the position vector of the forward trimmer panel 

! 

center of gravity is as follows: 


p « X” 1 " + Y" 1 " + z M 1 " , 

ccg eg x eg y eg z 5 


(ID 


where, because of symmetry, Y" = 0. This transfers to the forward trimmer 
hinge axis system (which is mutually parallel to the wing panel axis system) 
through the forward trimmer panel deflection angle, 6^. This can be transferred 
to the wing hinge axis system through the wing panel deflection angle, & . See 
figure 20. 

Before the forward trimmer center of gravity position vector can be 
written in the wing hinge axis system, the forward trimmer hinge axis (origin 
of both forward trimmer axis systems), must be located with respect to the wing 
hinge axis (origin of both wing axis systems), so that 

ccg^WH = p ccg + P CH = p xCH*xH + p zCH*zH ' 

The result is 

(p ccg>WH = [X "cg COs(5 c + 6 P> + Z ”ccg Sln(6 c + V + X ’CH COS 


+ Z' CH sin 6 p ]l vU + [-X"_ sin (6^ + <5 P ) - X' pH sin 6 P + Z ' rv cos (12) 


P J xH 


eg 


CH 


CH 


P J zH‘ 
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FIGURE 19. ILLUSTRATION OF AXIS SYSTEMS 


Canord center 
of gravity 


* CCG^ 


Wing /boom 




Aircrafl 
center of 
grovity 


FIGURE 20. ANGULAR RELATIONSHIPS OF AXIS SYSTEMS 


64 


Since the vertical offset of the forward trimmer center of gravity from its 
hinge is effectively zero, the second term in the preceding equation may be 
dropped. 

Now, define wing hinge velocity terms, V u and oo where V u is inertial 

H H ri 

■ velocity and u) is angular velocity of the wing hinge in that system, or 
H 

\ ' U1 XH + V V + > < 13 > 

and 

^ W H = P *xH + q ^yH + r ^zH ' 

Then the trimmer panel center of gravity inertial acceleration may be 
| found in the wing hinge axis system as follows: 

1 \cg = (P ccg ) WH + (15) 

: ^ccg = (p xCH + U) 1<H + V ^yH + (P zCH + W) *zH 


+ P xCH dt ( ^xH ) + P zCH dt ( ^zH ) ' 


*ccg = (P xCH + qP zCH + qP zCH + U) X xH + ( * P xCH + rp xCH 

PP zCH ~ P ^zCH + ^yH + ^ P zCH " qP xCH ~ qP xCH + ^zH 
+ ( ^xCH + qP zCH + U) dt ( ^xH ) + (rp xCH “ PP zCH +V) dt ( ^yH ) 


+ ( qP xCH + P zCH + W) dt ^zH^ ’ 
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Unit vectors differentiate to: 


dt (1 xH ) rl yH ” ql zH 


dt ^yH-* 


rl „ + pi „ 
xH zH 


dt ql xH pl yH * 

Therefore, the preceding expression for forward trimmer inertial acceleration 
becomes, upon expansion, 

<Vg>* - { -' 5 c r cg sin(4 c + 4 p> - «c X "cg cos(4 c + 4 p> - «cV" cos(4 c + V 
+ V^'cg sln<5 c + 4 p> - X 'cH 4 p + Z, CH 6 pl 

+ ‘pWc^cg COS<5 c + 4 p> " V*'cg cos < 4 c + S p) + X ' CH cos 4 p 

+ Z' CH Pin « P )] + (-r 2 - , 2 )[X" cg cos(6 c + « p ) + X' CH cop « p 
+ Z ’cH Sin 6 p ] + (q + pr)[ - r cg S±n (6 c + 6 p> " X 'cH Sin 6 p 


+ 2 1 cos <5 ] + 2 q [-6 X" cos (6 + 6 ) - 6 (X 1 cos (6 + 6 ) 

CH p ^ c eg c p 7 P eg c p 7 


c p' 


+ X 'CH COS 4 p + Z 'CH sl " 4 p» + " - rV + i*H 


‘Vg’, ■ t(l,< I + f)lx "cg COs(4 c + 4 p> + X ’CH COS S p + Z 'cH Sl " 4 p> 

+ 2r [-« c X "cg slnW c + 4 p> + 4 p <" *'cg sin(S c + 4 p> ‘ X 'CH aln 4 p 

+ Z ' CH cos 6p) ] + (-p + qr) [-X" sin (6 £ + 6 ) - X' CH sin 6 p 

+ Z' cos 6 ] - 2p £—6 X" cos(6 + 6 ) - 6 (x" cos(<5 + 6 ) 

CH p c eg c p p eg c p 

+ X' „ cos 6_ + Z' sin 6 ) ] + V + rU - pW} 1 „ . (19) 


(Equation 19 continued) 

pr) [X' 

ccg' z - ' eg 


= { (- 9 + P r M x "„„ cos ( 6 C + 6 p ) + X ’ CH cos 6 p + z ’ ch sin 6 p ] 


- 2q[- 6 X" sin (6 4 <5 ) 4 6 (- X” sin(<5 4 <5 ) 

^ c eg c p' p eg c p 

~ X, CH Sin 6 p + Z ’cH COS 6 p )] + (-q 2 -P 2 >[- X " C g sin ( 6 c + V 

- X’ CH sin 6 p + Z' CH cos 6 p ] - 6 c X" cg cos(« c + « p ) 

+ 6 Z X” sin (6 + 6 ) + 6 6 X" sin(<S + 6 ) - 6 [X" cos(6 

c eg c p' c p eg c p p eg c 


+ 6_) + X' CH cos 5 p + Z'o, ain S p ] - y-S c X" cg al»(« c + « p ) 


4 6 (- X" sin (6 4 6 ) - X’ sin 6 4* Z’ cos 6 )] 4 W 

p eg c p CH p CH p 


- qu + pv) i 2H • 


(19) 


Finally, collecting terms by sines and cosines yields: 


(a ccg ) x = ~ (<5 c + 6 p + q + P r )[ x " cg sin ( 6 c + V ] - [(6 c + V + 2q(6 c + V 


c P 


eg 


+ q 2 + r 2 ] [x"^ o cos(6 ( , + S n ) ] + (6 + q + pr)(-X' CH sin 6 


eg 


c P 


+ Z ’CH COS + K )2 + r2 ' (X 'cH C0S S o + Z ’CH Sl " V + » 


CH 


- rV 4 qW 


(a Ccg ) y = _[2r(<S c + V " P + qrltX "cg Sin(6 C + 6 p )] + [2P(6 C + 6 P ) 


+ pq + r] [X” cos(6 + 6 )] + (2r6 - p + qr) (-X' sin 6 


eg 


c P 


+ Z' CH cos 6 ) + (2p6 p + pq + r)(X' CH cos 6 p + Z ' C H sin V + V 
- pW + rU 


( 20 ) 
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(Equation 20 continued ) 

< a ccg> 2 * [<« c + «p> 2 + 2< 5 « c + « p ) + P 2 + <! 2 ]!X" I . g sln(S c + 6 p) ] 

- (6 + 6 + q - pr)[X" cos(6 + 6 )]-[(q + 6 2 ) -[(q+6 ) z 

u F c F p p 


+ P ^ ch sin 6 - + Z? ru cos - (<L + q “ P r)(X f ru cos 6 


CH 


CH 


+ Z’ „ sin 6 ) + W + pV - qU . 

Cn p 


( 20 ) 


Then, applying the fundamental Newtonian law, the three equations 
describing the forces at the origin of the wing hinge axis system that are 
associated with the acceleration of the forward trimmer panel are as follows: 

(F ) = m (a ) 

xH c c ccg x 

(F ) = m (a ) 

yH c c ccg y 


(F ) = m (a ) 

zH c c ccg y z 


( 21 ) 


Wing panel force equations . In the wing panel axis system in 
figure 19, the position vector of the wing panel plus boom center of gravity is: 


p Tm = X’ 1 f + Y ! 1 , + Z T 1 f . 

WBcg eg x eg y eg z 


( 22 ) 


Again, Y 1 Q g = 0 because the free body is symmetrical about the x’-z 1 
plane. The appropriate angular transfer must be made to write this position 
vector in the wing hinge axis system. 

(p Tn > ) Tru = (X 7 cos 6 + Z r sin 6)1+ (-X 7 sin 6 

WBcg WH eg p eg p xH eg p 


+ Z cos 6 )1 TT 
eg p zH 


(23) 


Simplifying momentarily to 


P WBcg P xH 1 xH + P zH 1 zH 


( 24 ) 


yields, with V and a) the same as before, 
rl H 
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(25) 


V, 


V WBcg * p WBcg^WH + V H’ 

(PvH + + + (P zH + + P^U 77 (Ku) + P„u 4^ 


WBcg H xH xH yH 


zH ' m xH dt v xH 7 T k zH dt v zH 7 


(26) 


^WBcg = (U + qp zH + qp zH + p xH U xH + (V " PP zH " PP zH + * p xH + ^xH^yH 


+ (W - qp xR - qP xR + P zH )l zH + (D + qP zR + P xH > IF (1 xH } 

+ *- V ~ PP zH + rp xH' ) dt ^yH-* + " qp xH + p zH^ dt * 


(27) 


Expanding this in its component form yields: 

( a_ __ ) = [-(q + 6 ) 2 - r 2 ](X' cos 6 + Z' sin 6 ) + (6 + q + pr) 

v WBcg 7 x ^ p eg p eg p p 

(-X' sin 6 + Z' cos 6 ) + U + qW - rV 

eg p eg p 7 

(a TT „ ) = (2p6 + pq + r) (X' cos 6 + Z' sin 6 ) + (2r6 - p + qr) (28) 

WBcg 7 y p eg p eg p p 


(-X' sin 6 + Z' cos 6 ) + V - pW + rU 

eg p eg p 


: s 2 


(a Tir) ) = (-6 + pr - q) (X 1 cos <5 + Z' sin 6 ) + [— (q + 6 ) 

WBcg 7 z p eg p eg p p 


- P ] (-X' sin 6 + Z' cos 6 ) + W + pV - qU 

F eg p eg p 

Again, applying the fundamental Newtonian law, the three equations 
describing the forces existing at the origin of the wing hinge axis system that 
are associated with acceleration of the wing panels are: 


^xH^WB = “wB^WBcg^x 
^ylPwB = m WB (a WBcg ) y 
( *zH*WB = “wB^WBcgh • 


(29) 


Trimmer- pan el moment equations . The trimmer- panel moment equations 
will follow the general form for moments about offset points: 

M = H + (p) x m(R) . (3°) 
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The first term will be referenced to the forward trimmer panel system 
where all inertias and angular velocities will be those of the forward trimmer 
panel about its hinge axis system. The position vector, p, will be expressed in 
the forward trimmer panel axis system initially, as will the linear momentum 

15 . 


term, mR. All of the terms will be transferred to the wing hinge axis system 
and combined at that time. 


momentum 

term is 

of the 

form: 

1 

XX c 

- I 

xyc 

-I 

xzc 

- V 

0) „ 

X 

, - I 

I xyc 

I 

yyc 

- I 

yzc 

: V 

; - I 

xzc 

i_ 

- i 

yzc 

I 

zzc 

! “ z "_ 


V 

1 y" 

V 


(31) 


Since the hinge axis system is composed of axes of mass symmetry, all products 
of inertia are zero. The angular velocity terms can be found by considering the 
relationship of the canard to the wing hinge axis. The only freedom the forward 
trimmer has with respect to this point is in pitch and is solely a function of 
the displacement angles 6 and 6 , such that 

p c 


co H 
xyz 


= V' 1 x" 


00 „1 „ 

y y 


„i ,■ 

Z Z 


(32) 


III) M 

xyz 


= p i u + (q + (S + 

xH n c 


6 > X u 
p yH 


+ rl 


xH 


(33) 


H q! can then be differentiated with respect to time to yield the first term in 
the moment equation. 

The position vector, which is the arm to compute the moment of 
momentum, runs from the forward trimmer hinge point to the forward trimmer 
center of gravity and can be expressed as: 


p o' 


= X" 1 „ + Z" 1 „ 
ccg eg x’ eg z 


To simplify this a bit, let Z" = 0 to reflect the nominal difference in the z" 

eg 

coordinate of the forward trimmer center of gravity. Transferring this to the 
wing hinge axis system through the forward trimmer panel and wing panel axis 
systems yields 


V- 


= X" cos (6 + 6 )1 „ - X" sin (6 

ccg eg c p' xH eg c 


+ 


6 )1 u 

p zH 


(34) 
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Next, the vector R , must be found. 
9 o 


R , * a + p f 
o’ o K o - o 

V = i + ^2 «VV + 2 ’o'V> 

at 


(35) 

(36) 


R , = V + — x [(X' , cos 6 + Z' , sin 6 )1 + (- X' , sin 6 

o o ^2 o p o p xH o p 


+ Z' , cos 6 )1 u ] 
o p zH 


(37) 


Now, for convenience, let 


p , = X' , cos 6 + Z' , sin 6 

ox o p o p 

p , = - X' , sin 6 + Z' , cos 6 

o z o p o p 


(38) 


Then the vector becomes: 


R , = (U + p , )1 _ + VI _ + (W + p , )1 „ + (U + 2p , ) 4~ (1 u ) 

o' ' K o x xH yH o'z zH v o'x dt xH 

2 2 

+ V — (1 _) + (W + 2p , ) — (1 „) + p , (1 H ) +p , (1 „) . 

dt yH o z dt zH o x xH 7 o z j4 .2 zH 


dt 


dt 


Expanding yields, 

R , = {-[(6 + q) 2 + r 2 ](X' , cos 6 + Z' , sin 6 ) + (6 + pr + q) 

ox p opopp 

(- X' , sin 6 + Z' , cos 6 ) + D + qW - rV} 1 _ 

o p o p xH 


(39) 


R , = [(2p6 + pq + r) (X' , cos 6 + Z' , sin 6 ) + (2rS + qr - p) 

o y p o p o p P 


(40) 


(- X’, sin 6 + Z' , cos 6 ) + V + rU - pW]l „ 

o' p o p yH 


it 

R , 


O Z 


{-(6 - pr + q) (X f , cos 6 + Z T f sin 6 ) - [ (5 + q)^ + p^] 

p r o P° P P 


(- X 1 , sin 6 + Z* , cos 6 ) + W 4- pV - qU} 1 


zH * 
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Finally, since the mass involved is that of the forward trimmer alone, the 
moment equations become, 


(M f ) = { - p(8 + 8 + q) I + rl - m X 1 [ [- X , (2ro + qr - 

o’ X L v v c p H XX zz c cg LL o’ p M 

^ c c r 

+ Z 1 , ( 2p<5 + pq + r) ]sin 6 + [X 1 f (2p8 4- pq + r) + Z’ f (2r8 

o P p o p o p 


4* qr - p) ] cos 8 + V + rU - pW] sin (6 4-8)4- [pi 4- r(6 

P J P XX c 

+ 8 + q) I ] cos(5 4- 6 ) - r(8 + 6 4- q)I }1 

p zz^ c p c p yy c xH 

(M.) = ( \t> 2 I + r 2 I + m X" [[X' ,(6 + q + pr) + Z' ,[(6 

o y \ [_ xx^ zz^ c eg o p o p 


O’ p 


+ q) 2 + r 2 ]] sin 6 + [X' , [(6 + q) 2 + r 2 ] - Z’ , (6 + q 

p o p n o p 

+ pr)] cos <5 - U + rV - qW]] sin(6 + 6 ) + [pr(I -I ) 

p J c p [/ xx zz 

c c 

+ m c X "cgih x, 0 »t(6p + q) 2 + p 2 ] + Z' o ,(6 p + q - pr)] sin 6^ 

+ [X' o ,(6 + q - pr) + Z* o , [(6^ + q) 2 + p 2 ]]cos + W - qU 


+ pV] cos (6 + 6 ) + (6 +6 + q)I }l TT 

J C p C p yy c yH 


(M t ) - { [- pi - r(6 + 8 4- q)I ]sin(8 4-6)4* j— p(6 4- 8 

o' Z XX c c p ^ ZZ C C P L c P 


+ q) I + rl + m X n [[- X’ f (2r8 - p + qr) 4- Z’ f (2p6 

xx zz c ce u o p o r p 

c c e 

+ r 4- pq)]sin 8^ 4* [X’ ol (2p8^ 4* f 4- pq) 4- Z’^, (2r8^ ~ P 


4* qr)]cos 6 + V + rU - pW] cos(8 4- 8 ) 

P J op 


4- p (6 4 8 + q) I }1 



! Wing/ boom moment equations . The wing/boom moment equations are 

written most conveniently in the panel axis system (shown in figure 19), 
because in this system, the moments and products of inertia are constants. The 
! moments are then transformed to the hinge axis system for later use. 

An unusual feature of the panel axes is that the origin is displaced 
t from the wing/boom center of gravity. Because of this, the more general form 
of the principle of the consdervation of moment of momentum must be used. This 
is: 

-x (42 

M = H + (p) x m(R) 


The components of the H vector are the inertial terms found in the 
t conventional Euler equations for the rotation of a rigid body. These are not 
I rederived here because they are developed in many texts. Their form is more 
general than that shown for the forward trimmer in equation (31) in that 

! 

products of inertia are included. 

| The second term, caused by the offset center of gravity, requires the 

development outlined below. 

In the wing panel axis system, the position vector to the wing/boom 
center of gravity, p, is constant, and is given by 


p o-WBcg X WB 1 X' Z WB 1 Z' 

I 

In the wing hinge system, this is 


(43) 


p = (X' cos 6 + Z" sin 6 )1 + (-X' sin 6 + Z' cosS )1 

O-WBcg WB p WB p xH WB p WB 


p' zH 

(44) 


The angular momentum term will be, then, the time derivative of 


H 0= 


\i 


^xxWB byWB bzWB 


byWB I yyWB ^yzWB 


-t - I I ! 

xzWB yzWB zzWB J 


xH 


U) 


yH 


zH 


.i; 


Finally, the linear momentum can be found by noting that 

R = U1 „ + VI „ + W1 


xH 


yH 


! 1 zH 


xH 


yH 


zH ' 


(45) 
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R o = (U + qW - rV)l xR + (V + rU - pW)l yR + (W + pV - qU)!^ . (4f 

The moment equations can be written in their final forms by crossing 

-V *1? 

3 0-WBccg int ° V multipl y in g b y “yg’ and adding H q to obtain: 

<Vx = ^ I xxWB + (pr " q ' )I xyWB ” ( ' Pq + r ^ I xzWB “ qrI yyWB + ^ q ) I yzWB 
+ qrI zzWB " “WB + rU - PW)( " X ’ W B Sin 6 p + Z ’WB C ° S 6 p )]t xH 
<Vy = [prI xxWB " (qr + p)I xyWB + (p “ r ^xzWB + qI yyWB + (pq ” r)I yzWB 
prI zzWB~ "Vb + pV ” qU ^ X 'wB C ° S 6 p + Z WB Sln 6 p^ 

(47 

+ (U + qW - rV)(-X’ WB sin 6 p + Z’ Wfi cos 6 p ) ] i yR 

<Vz = PqI xxWB + (^ " p ^ I xyWB + < ‘ qr ” p)I xzWB + PqI yyWB 

' (pr + q)I yzWB + H zzWB + “WB + rU “ pW)(X ' W B COS 6 p 


+ Z ’WB S±n V ]1 zH • 

Fuselage moment equations . Consider the fuselage as a weight sus- 
pended from a wing and connected to it by a single pin joint so that it is free 
to rotate about this point unless constrained by physical stops or aerodynamic 
forces on some auxiliary trimming surface attached to it. It should be 
apparent, then, that longitudinal motion will occur about this wing hinge axis 
for the fuselage as well as for the free wing/free trimmer combination. To 
describe the moments acting on the fuselage, it is again necessary to consider 
the general form of the equation for the conservation of angular momentum. 

M p = H + (p f ) X m F (R Q ) . (48) 

Recall that there are three axis systems fixed at the intersection 
of the wing hinge line and the aircraft plane of symmetry: the wing panel axis 

system fixed in the wing, the wing hinge axis system mutually parallel to the 
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aircraft stability axis system, and the fuselage axis system with Xp positive 

forward and parallel to the fuselage reference line and Yp positive right and 

coincident with the wing hinge line. Z„ is orthogonal to these two axes, 

r 

positive downward. Both the fuselage center of gravity and the aircraft center 
of gravity are fixed points in the fuselage at any point in time and, hence, 
have fixed coordinates in the fuselage axis system. Assuming that the fuselage 
center of gravity lies in the aircraft plane of symmetry, the position vector 
from the wing hinge point to the fuselage center of gravity will be: 


P F X Fcg X xF + Z Fcg 1 zF 


(49) 


Since the linear velocity of the origin of the stability axis system is the 
velocity of the fuselage center of gravity, 

K * ut x + vt Y + “z ' 


(50) 


However, the wing hinge axis system is parallel to the aircraft stability axis 
system, so 


and 


R o - U1 XH + V1 yH + W1 ZH • 


K = ^xH + ^yH + ^zH + * <V • 


(51) 


(52) 


where 


W H = pl xH + ql yH + rl zH * 


Then, the second term in equation (48) becomes: 


(p F )x m£R) = - m F Z Fcg (V - P W + rU)l xH + m f [Z fcg (U + qW - rV) 


(53) 


- X fcg (W + P V - qU)]l yR 4- m F X Fcg (V - P W + rU)l zH . 

The remaining term on the right side of equation (48) is, as before, the 
inertial term found in the conventional Euler equations for the rotation of a 
rigid body. 
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become: 


The components of the applied fuselage moment defined by equation (48] 


M xF " ^xxF + ( J> r - ^ + + (f 2 - q 2 )I yzp 

+ ,r(I zzF - I yyF ) - VFcg^ - P" + r«) 

V ’ Pl(I xxF - r z*F> - « r + » V + (p2 - r2>I xzF + P P yy F 

yy (54) 

+ (pq - r)I yzF + m F [-X Fcg (W + pV -qU) + Z pcg (U + qW - rV) ] 

M zF ' P ’ (I yy F - W + - P 2 )V + (IF - P)I XZF j 

- (Pr + q)I yzF + rl zzp + n^X^V - P W + rU) I 

! 

The moments applied to the fuselage assembly, represented by the sides j 
to the left in equations (54) , contain contributions from the reversed effective 
forces and moments of the wing panels, booms, and trimmer. In actuality, they 
also contain gravity moments due to the weight of the fuselage, wing panels, 
booms, and trimmer* 


Summation of Equations 

Force equations . Summation of the force equations for the total air- 
craft will be simple and of the form: 


’ + F , . - m g sin 6 = m (U + qW - rV) 

x, aero x, thrust a fo a v n 

+ F . , . + m g cos 0 sin <b = m (V + rU - pW) 

y , aero y, thrust a B y a v v ' 


(55) 


F + F ,_ + m g cos 0 cos 6 = m (W + pV - qU) 

z , aero z, thrust a & Y a v v n 

Moment equations . The moment equations will express that the sum of 
all externally applied moments, gravity moments, and moments generated by the 
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wing /boom/ trimmer combinations on the fuselage must equal the fuselage inertial 
moments, or 


Z M , . , . = H _ + (p ) x m_(R_) 

applied to F 1 F F F 

fuselage 

Broken into components, the scalar equations are: 


( 56 ) 


M x,aero + M x, thrust + M x, fuselage + M x, wing/boom/ trimmer 

gravity gravity 


JL M 4* M 

x,due to reversed effective x, reversed effective forces 


moments at the hinge 


(inertial) at the hinge 


( 57 ) 


= pI xxF + (pr - q)I xyF - (pq + r)!^ + (r 2 - q 2 )I yzF 


+ qr(I _ - I _) - m (V - pW + rU) . 

M zzF yyF a Fcg ^ 


m +M +M +M 

y,aero y, thrust y, fuselage y, wing /boom/ trimmer 


gravity 


gravity 


+ ^y, reversed effective + ^y, reversed effective 


moments at hinge 
"yyF ' _/I yzF 


= ql,„. w + (pq - - (qr + p)I to17 + (p 2 ~ r 2 )I xzF 


forces at hinge 
r xyF 


( 58 ) 


+ pr(I xxF - I zzF ) - m a t* Fcg (W + pV ' qU) - Z Fcg (il + qW - rV) ] . 


z,aero + z, thrust z, fuselage z , wing /boom/ trimmer 

gravity gravity 


+ M + m 

z, reversed effective z, reversed effective 


moments at hinge 


forces at hinge 


( 59 ) 


= tI zzF + (qr ' p)I xzF “ (pr + q)I yzF + (q2 ' p2)I 


xyF 


+ pq(I - I ) + m X (V - pW + rU) 
yyF xxF a Teg 
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M and M - terms are dealt with in appendix B, and the gravity force 

aero thrust & 

terms appear below: 


x, fuselage 
gravity 


'Fcg m F 8 C0S 


0 sin <J> 


M 

y, fuselage 
gravity 

M z, fuselage 
gravity 


■ C0S d cos <t> - Z Fcg m F § Sin 0 



( 60 ) 


In deriving the reversed effective moment terms in equations (57) to 
(59), note that the y-components of these moments cannot be transmitted to the 
fuselage through the wing hinge, since by definition, the wing hinge is a 
frictionless pin joint. Practically speaking, friction forces would exist, but 
the calculation of their magnitude is beyond the scope of this work. The 
reversed effective moments at the hinge, then, will be: 


Mx, reversed effective ^ox + ^o T x^ 

moments at hinge 

m , -- = - (M + M , ) 

reversed effective oz o z 

moments at hinge 
The total equation is: 


M = - (M + M t )1 u + 
ox o T x xH 


01 yH- < M o. +M o'* ) 1 zH 


(61) 


(62) 


Next, consider the moments due to reversed effective forces at the 
hinge. The inertial accelerations of the trimmer center of gravity and of the 
wing/boom center of gravity were expressed in the wing hinge coordinate system. 
The resulting reversed effective forces are, therefore, expressed in that 
system. These forces, in fact, act at the center of gravity of these two free 
bodies. It is possible to show how, with only one exception, they can be 
accounted for as acting at the hinge points of the trimmer and of the wing. 
Consider the explanation, by way of figure 21, as follows: 
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FIGURE 21. FORCES AND ACCELERATION VECTORS ACTING ON 
WING /BOOM/ FORWARD TRIMMER ASSEMBLY 


To the extent that these derivations have considered the "reversed 
effective moment" due to center of gravity offset in the moment-derivations 

-v IS. 

(p) x m(R) about the hinge points, it would be redundant to consider it again. 
Therefore, both of the reversed effective forces can be shifted to their 
respective hinge points. For the wing/boom force, this is a complete step — the 
reaction is acting at the hinge and the moments resulting on the fuselage are 
calculable. For the trimmer force, the process is more involved. First, the 
force can indeed be viewed as acting at the trimmer hinge. Furthermore, it 
can then be shifted to the wing hinge as long as the gross moment due to this 
shift is also accounted for. 

The process of getting the reactions at the wing hinge is, therefore, 
straightforward, except for this one moment. Note here that the y-components 
of this moment, as for the pure inertial moments, cannot be transmitted through 
to the fuselage. 

In summary, the reversed effective force acting at the wing hinge is 

given by: 


^hinge ^ F x, trimmer , H + F xH^xH ^y, trimmer, H + F yH^yH 


- (F . + F ) 1 

z, trimmer, H zH zH 


(63) 
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The reversed effective moments due to forces at the hinge, then, will be this 
force multiplied by its moment arm. In vector notation, this is: 


M . = (p) x (F ) > 

reversed effective hinge 

forces at hinge 

where p~ is the position vector of the wing hinge axis with respect to the 
summation point — the wing /boom center of gravity: 


( 64 ) 


» - X1 x H + X1 zH ' 0 ' 


(65) 


Then, the moment due to reversed effective forces at the wing hinge becomes: 

(p)x(p, . )= Z(F „ . „ + F U )1 u + [X (F . . u + F „) 

hinge y, trimmer, H yH xH z, trimmer, H zH 


- Z(F „ . u + F u )]l u - X(F _ . u 

x, trimmer, H xH yH y, trimmer, H 


( 66 ) 


+ ^H^zH 


0 


The components of this yield partially the M reversed effective forces for 

x , y , z 

equations (57), (58), and (59). What is now necessary is the contribution from 
the cross-product of the position vector from the wing hinge to the trimmer 
hinge into the reversed effective force acting at the trimmer hinge. The force 
is just: 


p = _ p 1— F 1 

trimmer hinge x, trimmer, H xH y, trimmer, H yH 

reversed effective 

force (67) 


F 1 

z , trimmer ,H zH 


The necessary position vector is constant in the wing panel axis system, but it 
must be expressed in the wing hinge axis system. In the panel axis system, the 
vector is: 


P CH " X ' CH 1 x' 


+ Z ' CH 1 z' 


( 68 ) 


The final form of these moments, due to the total reversed effective forces, is: 


80 



M [Z(F . . „ + F ) + (Z'„ cos 6 

x, reversed effective y, trimmer, H yH CH p 

forces at hinge 

- X f sin 6 )F . 1 

CH p y, trimmer, H xH 

M = f x(F + F ) __ z (f 

y, reversed effective z, trimmer, H zH' ' x, trimmer, H 


forces at hinge 


+ f kh>1V 


M , .. = [- X(F „ . + F ) - (X' „ cos 6 

z, reversed effective y, trimmer, H yH CH p 

forces at hing 

+ Z' „ sin 6 )F U ]1 „ . 

CH p y,tnmmer,H zH 


( 69 ) 


Last, the gravity moment must be calculated. The components of the 
gravity force are shown in equations (55). Defining a position vector from the 
wing hinge to the wing/boom/trimmer center of gravity as follows: 


ms = QC to cos 6 + Z s in 5 )1 „ + (-X 

-WBccg WBccg p WBccg p xH Vi 

+ Z~ _ cos 6_)1_„ , 


the moment due to gravity at the wing hinge is, with 


W, 


WBC 


= -m WBC 8 sin 6 hH + "W 8 cos 0 sin * 1 


yH 


+ m WBC g cos 0 cos <j> 1 


zH 


(p) x '"WBC' " ‘ (£ WBC COS S p - *WBC Sl ” S p )m HBC 8 C0S 9 Sin * 1 


xH 


(71) 


+ (non-transmittable y-terms) 

+ ^¥BC COS 6 p + W Sin 6 p )m WBC 8 C ° S 9 Sin + hH • 

Finally, all these parts may be summed to yield the total moments 
acting on the fuselage at the wing hinge due to inertial forces: 
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M 


x,aero 


+ M X, thrust ^Fcg 11 ^ + ^WBC C ° S 6 p 


^BC Sin 6 p )m WBC ]g COS 6 Sln * 


(M + M , ) + [Z (F „ . + F ) + ( Z' cos 6 

ox o'x / L v y, trimmer ,H yH CH p 


X CH sin ‘ 5 p ) F yjtrlmmer jR ] P I xx p + (P r ^^xyF ( ' Pq + ^xzF 
+ ^ - 0 2)I yz F + ^zzF - VF* - m a Z Fcg^ ’ t>« + rD) ‘ 


M 


+ M 4- “ X-, m^g cos 0 cos <J>- Z m^p- sin 0 + [X(F . TT 

y,aero y, thrust teg F e Y Fcg F^ z, trimmer, H 


+ F ) - Z(F . + F ) ] = ql _ + (pq - r)I ^ - (qr + p)I _ (72) 

zH x, trimmer, H xH yyF yzF M v xyF 


+ ( P - r )I xzF + prd xxF - I zzp ) - m a [XFcg (W + P V - qU) - Z^ (U + qW - rV) ] 


Fcg 


M 


z ,aero 


+ M z , thrust + [X Fcg m F + <*WBC C0S S p + Z WBC Sin V"wBC lg C0S 9 sin 


- (M + M , ) + [- X(F _ . „ + F TT ) - (X’ cos 6 - Z’ „ sin 6 ) 

oz o'z y, trimmer, H yH CH p CH p' 


F ] = rl „ + (qr - p)I „ - (pr + q)I + (q^ - p^)I „ 

y, trimmer, H zzF M v xzF VF M yzF H v xyF 


+ pq (I -I -) + m X-. (V - pW + rU) 
yyF xxF a Fcg v 


M , M , M are given in their entirety in equations (47); M , , M , , M , are 
ox oy oz 6 J M * o x o y o z 

given in equations (41); trimmer forces are expressed in the wing hinge 
coordinate system in equations (21); and, finally, wing/boom force equations 
are given in the wing hinge system in equations (29). 
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APPENDIX B. DEVELOPMENT OF LINEARIZED 
LONGITUDINAL EQUATIONS OF MOTION 


Introduction 

The system is composed of three rigid bodies coupled together by 
geometric and kinematic constraints. The three bodies are the free trimmer 
panel, the wing/boom assembly, and the fuselage assembly. 

Force and moment equations are developed for each free body in 
linearized form, with motion confined to the vertical plane. Although 
these individual equations implicitly contain the geometric and kinematic 
constraints, the forces transmitted through the trimmer hinge and the wing 
hinge are treated as separate, dependent variables. Furthermore, because of 
the complexity of the representation of unsteady aerodynamic forces, it is 
also convenient to treat the circulatory lift coefficients of the wing and 
trimmer surfaces as additional explicit variables. 

As derived and employed in this study, the equations represent a 
departure from conventional analysis. Although the number of dependent equa- 
tions could have been reduced, and all of the forces and moments could have 
been referred to the total aircraft center of gravity, as is customary, the 
physical significance and origins of the individual terms in the equations 
would have become obscure. In short, the individual equations are retained 
as they were derived to facilitate their verification — not a trivial con- 
sideration for such a complex system. 

Coordinate Systems 

Four coordinate systems were employed in the derivation of the linear 
longitudinal equations. These are identical to the systems used in the non- 
linear derivation, but they are briefly described here for completeness. 

Wing hinge axis system . The wing hinge axis system, x^, y^, z^, 
is the primary axis system used to define the fuselage pitch angle, pitching 
rate, and translational velocity components of the total aircraft motion. 

The origin is at the intersection of the wing hinge axis and the plane of 
lateral symmetry of the aircraft. The x^-axis is positive forward in the 
direction of flight, in the equilibrium condition. The y^-axis is positive 
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toward the right wing. The z^-axis completes the right-handed set, positive 
downward. 

The wing hinge axes are fixed with respect to the fuselage assembly 
and rotate with the fuselage in pitch. 

Wing panel axis system . In the wing panel axis system, x f , y f , z f , 
the origin is coincident with the origin of the wing hinge axis system, but 
the panel axes are fixed to the wing and rotate in pitch with the wing panels 
The x ! -axis is parallel to a reference chord line in the wing and is positive 
toward the leading edge. 

Trimmer hinge axis system . The trimmer hinge axis system, x^, 
y^, z^, is parallel to the wing panel axis system, but the origin is at the 
intersection of the free-trimmer hinge axis and the plane of lateral symmetry 

Trimmer panel axis system . The trimmer panel axis system, x", y", 
is coincident with the origin of the trimmer hinge axis system, but is fixed 
in the free-trimmer surface and rotates in pitch with it. The x"-axis is 
parallel to a reference chord line in the trimmer and is positive toward the 
leading edge. 

Symbols 

acceleration of center of gravity of free trimmer, 

2 2 
meters/second (feet/second ) 

acceleration of fuselage center of gravity, meters/ 

2 2 
second (feet/second ) 

acceleration of the center of gravity of wing boom 

2 2 

assembly, meters/second (feet/second ) 
mean aerodynamic wing chord, meters (feet) 
mean free trimmer chord, meters (feet) 

fuselage drag coefficient 


a 

ccg 

Fcg 

a WBcg 

C 

C 

c 



C T , C. 


profile drag coefficients of free trimmer and wing, 
respectively 

lift coefficients of free trimmer and wing, respectively 




w 




t 


equilibrium lift coefficients of free trimmer and wing, 
respectively 

C lift-curve slopes per radian for free trimmer, wing, 

L 

Ct 

F and fuselage, respectively, per radian 
= 9 c /9<5 t , per radian 


C 

m 

C 



C 


m 


e 


C 


m 


q 


c 



pitching moment coefficient 

= 9c /9a_, per radian 
m F 

= 9C /9e per radian 
m 

= 9C / 9 ! -5“ | , per radian 
m • 2U I 

’ o } 

= 9C /d& , per radian 
m t 


dC D 



E 

F , F 

x , z , 
c~>b c-*b 

F , F 

Vhf z w->-f 


g 



i 

c 


induced drag factor 

ratio of wing semiperimeter to span 

components of force transmitted from free trimmer to wing 
boom, newtons (pounds) 

components of force transmitted from wing to fuselage 
newtons (pounds) 

2 2 

acceleration of gravity, meters/second (feet/second ) 

transfer functions relating lift coefficient to angle of 

attack for free trimmer and wing, respectively 

transfer functions relating lift coefficient to gust 

velocity for free trimmer and wing, respectively 

pitching moment of inertia of fuselage assembly about wing 

2 2 

hinge axis, kilogram-meters (slug-feet ) 

pitching moment of inertia of wing/ boom assembly about wing 

2 2 

hinge axis, kilogram-meters (slug-feet ) 

pitching moment of inertia of free trimmer about its hinge 

2 2 
axis, kilogram-meters (slug-feet ) 

distance from quarter chord of wing to quarter chord of 
free trimmer, measured positive forward, meters (feet) 
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I 

t 


V "W “f 


q 

s 

s 

c 

U 

u 

V 

g 

W 


X 

ac 


X F 

A A 

X , X 
c 



X" 

eg 


Z 


F 


7 ' 

CH 


Z ' 

eg 


V W 


6 

c 


distance from quarter chord of wing to quarter chord of 
fuselage mounted horizontal tail, measured positive 
rearward, meters (feet) 

mass of free trimmer, wing/boom assembly, and fuselage 

assembly, respectively, kilograms (slugs) 

pitch rate of fuselage, radians /second 

2 2 

wing reference area, meters (feet ) 

2 2 

free trimmer reference area, meters (feet ) 

true airspeed, meters/second (feet/second) 

increment in airspeed divided by equilibrium airspeed 

gust velocity, positive upward, meters/speed (feet/second) 

component of velocity along z^-axis, meters/second 
(feet/second) 

distance from leading edge to aerodynamic center of wing, 
meters (feet) 

distance from wing hinge forward to fuselage center of 
gravity, meters (feet) 

distance from hinge point to one-half chord point for free 
trimmer and wing, respectively, meters (*°ct) 
x ? -coordinate of trimmer hinge axis relative to wing hinge, 
meters (feet) 

x’ -coordinate of center of gravity of wing/boom, relative 
to wing hinge, meters (feet) 

distance of free-trimmer center of gravity forward of hinge 
line, meters (feet) 

z^-component of displacement of fuselage center of gravity 
with respect to wing hinge, meters (feet) 

z T -coordinate of trimmer hinge axis relative to wing hinge, 
meters (feet) 

z T -coordinate of center of gravity of wing/boom, relative 
to wing hinge, meters (feet) 

angle of attack of free-trimmer, fuselage, and wing, 
respectively , radians 

displacement of free“trimmer with respect to wing, positive 
leading edge up, radians 
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! 

■ e 
0 


X 

P 


P x » 


eg eg 


p 


*H Z H 
c c 


displacement of wing with respect to fuselage, positive 
leading edge up, radians 

displacement of control tab on free trimmer, positive 
trailing edge down, radians 
downwash angle, radians 

pitch angle of longitudinal fuselage axis with respect to 

horizon, radians 

Laplace operator, per second 

3 3 

atmospheric density, kilograms /meter (slugs/foot ) 
coordinates of wing boom center of gravity in wing hinge 
axis system, meters 
coordinates defined by equation (97) 


Trimmer Free Body 


The complete nonlinear expression for the acceleration of the center 
of gravity of the trimmer, expressed in the wing hinge axis system, is 
developed in appendix A. Eliminating the lateral-directional variables, the 
nonlinear expressions for the longitudinal acceleration are 

(a ) = -(6 +6 + q)[X" sin (6 + 6 )] - [ (6 + 6 ) 2 + 2q (6 + <5 ) + q 2 ] 

ccg x c p ^ eg c p c p c p 

fX" cos (6 + 6 )] + (6 + q) (-XL, sin 6 + Z' cos 6 ) 

L eg c p J p n CH p CH p 


- [(q + 5 ) K X q H cos <$ + Z^ H sin 5 p ) + U + qW 

(a ) = [(6 + 6 ) 2 + 2q(6 + 6 ) + q 2 ][X" sin (6 + 5 )] - 
ccg z c p M c p H eg c p 


(73) 


(6 c + S p + q)[X^ g cos (« c + S p )] - (q + i p > [-X' H sin S p 
+ Z €H COS S p ] ” (5 p + [X CH COS S p + Z CT 5in V + “ ‘ < i u 


In the equilibrium state, all angular rates and accelerations are 
zero, as are all translational accelerations. Furthermore, the x-axis of 
the wing hinge system is aligned with the direction of motion* 

Taking the differentials of equations (73), and denoting the 
equilibrium values of the variables by the subscript zero, the linearized 
form of the expressions for the acceleration of the trimmer center of gravity is 
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(a ccg ) x = U + [Z CH COS 6 p o ■ X CH Sin 6 p o ](4 + 6 p> 

+ [-X^ sin(6 + 6 c )](q + 6 + 5^) 

*0 0 V - 

< “ccg ) z ' " + '-"o 1 ' 1 + Sl ” 4 p - X CH “ S S p 1 * 

o r o 


(q 4- 6 ) 4- [“X cos (5 + 6 ) ] (q 4* 6 + 6) 

P c g P Q c 0 P c 


In equation (74) and in all subsequent linear equations, a variable 
without the zero subscript is understood to be a perturbation from the 
equilibrium value. 

The forces acting on the trimmer surface to cause the acceleration 
defined by equation (74) are composed of aerodynamic forces, weight forces, 
and the force transmitted to the trimmer through its hinge. 

The basic force equation for the trimmer is, then, 

F + m g 4- F, = m a (7f 

aero^ c b->c c ccg 

The circulation lift on the surface is based upon the angle of 
attack as defined at the one-half chord point by the velocity component nor- 
mal to the surface. 

In the trimmer panel axis system, 

V 

tan a = - — (7C 

c V „ 
x 

where V „ and V „ are the components of the inertial velocity at the one-half 

z X 

chord point. 

Differentiating equation (76), 


d(tan a ) = [dV „ - (tan a )dV „] 

C V f | z ox 

X 

o 


V „ = U cos (6 + 6 ) 

x o pc 

o r o o 


tan a 


sm (6 + 6 ) 

P c 

12 o_ 

c cos (6 4-6 ) 

o pc 

o o 
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i 


Substituting equation (78) into (77), 


dV 


d(tan a ) = 


sin (6 +6 ) 

q__ 

cos (6 + 6 ) 

P c 
o o 


dV 


x 


U cos (6 4- 6 ) 

o pc 

o o 


(79) 


The velocity components at the one-half chord point can be shown 


be: 

V x „ = [U + (q + 6 p )(-X^ H Sin 6 p + cos S p ) ] cos(6 p + S q ) 

+ [W - (q + « p )( X c H cos 6 p + Z ch Sln l5 p )][- sin ( 6 p + 6 C U 
[ V z „ = [U + (q + 6 ) (-XJh sin 6 p + cos 6 p )]sin(6 p + 6 c ) (80) 

| + [W - (q + <5 p M X c H cos 6 p + Z CH Sin 6 p Ucos(6 p + 8J 

^ -X (q + 8 + 8 ) 

i C p C 

f 

Taking differentials, 

dV „ = -U sin (6 + <5 ) (6 + 6 ) + cos(S + 6 ) [dU 

i x ° P G C o P c P o C o 

1 + (-X' sin 8 + Z' cos 6 ) (q + 6 )] - sin(6 + 8 ) x 

j CH P o CH P o P P o o 

! [dW - (X- cos <5 + Z' sin « )(q + « >] 

*0 o 

! 

i dV „ = U cos (6 +6 ) (6 + 6 ) + sin(6 + 6 ) [du (81) 

[ z ° P 0 C o P C P o C o 

' + ( -X c H sin 6 p + Z ch C ° S 6 p > ( q + l $p)]- cos ( <5 p + 6 C ) x 

[ dW - (X’ H cos 6 + Z’ H sin 6 ) (q + 6 )] 

r o o 

i - 

- X (q + 6 + 8 ) 

c p c 

1 , dW 

Substituting equation (81) into (79), and recognizing that — = « , 

o 

ihe angle of attack of the trimmer surface with respect to still air is 
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a = 5 + 6 + a_ - 

c p c f 


X' cos 6 + Z' sin 6 + X cos(<5 + 6 ) 

CH p CH p c p c 

]_o ]_o o 

U 


r X 


da 


( 82 ) 


(q + 6 ) + | - — cos (6 + 6 ) 6 + 

^ "n U ‘ n r» • i r* 


c [ c I dC L 
o J L W 

U W J 


The last term in equation (82) is the induced angle of attack at the 
trimmer caused by the vortex system about the wing surface. It is evaluated 
using a vortex-lattice computer program for the particular wing-trimmer geome- 
try being considered. 

Note that no explicitly defined lag is included to allow for the in- 
duced velocity to propagate upstream (fwd trimmer) or downstream (aft trimmer) 
following a change in the bound vortex system on the wing. On the other hand, 
the C itself lags the change in angle of attack by a time constant depen- 
dent upon the ratio of wing chord to flight speed. Consequently, the induced 
angle of attack increment at the trimmer will lag any change in the wing 
angle of attack by the same time constant. 

Equation (82) may be written 


where 


a = a + 6 +6 + [Cl]6 + [C2]q + [C3]6 + [C4 ]C T 

c f p c p LJ C 


X* cos 6 + Z’ sin 6 + X cos (6 +6 ) 


(83) 


CH 


CH 


Cl = C2 = — 

X 

C3 = - yp- cos(6 +6 ) 

U pc 

o o o 


U 


(84) 


da 


C4 = 


dC 


W 


In operational notation, where q = 0 and X is the differential (LaPlace) 
operator, 

[— l]a + [l]a + [1 + (Cl) X] 5 + [1 + (C3)A]S 4- [(C2)A]8 + [C4]C = 0 

c 1 p c L W 


(85) 


Equation (85) is the first in the set of linear equations describing 
the total system. It represents the angle of attack increment of the free 
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trimmer caused by motion of the aircraft with respect to still air. The incre- 

| ment in angle of attack caused by a vertical gust is not included in equation 

i 

(85) because the transient aerodynamic phenomena for gust penetration are 
different than for aircraft motion with respect to still air. For convenience, 
the vertical gust effects will be incorporated directly into the equation for 
| trimmer lift coefficient. 

The circulation lift on the trimmer surface arising from aircraft 
motion is based on the angle of attack of equation (85) plus an incremental 
camber caused by pitching and an effective increment caused by control tab 
deflection. The transient aerodynamic calculations and procedures are based on 
> reference 5, and a discussion of the techniques as applied to the previous 
^ analyses is contained in appendix B of reference 1. 

The circulatory lift caused by aircraft motion is 


C L = G 1 (t) [a c + ^T“ (8 p + q + V + C~ 6 t ] 

C _ c o r L 

mot. a 


( 86 ) 


G-. (t) is the lag function approximating the delay in the buildup 
of circulatory lift following a change in angle of attack. From reference 5, 
the lift function for an aspect ratio of 6, following a step change in angle 
of attack, is 

U 

-0.598 t 

C T = C. [1 - 0 . 361e ]Act (87) 

L . ij 

circ. a 


The corresponding transfer function relating lift coefficient to 
angle of attack is obtained by taking the LaPlace transform of the time 
derivative of equation (87) (ref. 6). 

The desired transfer function is 

G, (X) = 1 °- 361 - L c T (88) 

C L X + 0.598 “c 

c 

The circulation lift on the trimmer arising from a vertical gust 


is 


V 

C L = G 2 < t) u £ (89) 

C CO 

gust 
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G 2 (t) is the function approximating the lag in circulatory lift 
build-up foliowing the penetration of a vertical gust. Using the approxi- 
mate indicial response function from reference 5, and following the procedure 
described above, the corresponding transfer function, in Laplace operator 
notation, is 


G 2 (X) = C L [1 


0.4881 

U 

1 + 0.455 =2- 
c 


0.272X 0.193A ] 

U U J 

X + 1.04 jr- X + 4.71 =5- 

c 


(90) 


Since it is desirable to define the instantaneous gust velocity at 
the main lifting surface, the difference in gust velocity between the trimmer 
and the wing will be approximated by 



V i 

U U 
o o 


(91) 


Combining the various contributions the circulatory lift coefficient 
of the trimmer is 

C T 


C L = G (t) [a c + -- (6 + q + 6 ) + ^ + G 2 (t) ^ 

cc o r L. cooo 


(92) 


a 


Letting 



(93 


and multiplying by the denominator of equation (88), the final expression, 
in operational form is: 



u u 

[-(A + 0.598 ■=— • )]C + [0.639C t A + 0.598 C T ]a 

L L r. L C L C 

c c a c a 

c c 

7 U n 9 

+ [C L (C5)(0.639A Z + 0.598 -g— )]<5 +[C L (C5)(0.639A / 


U _ U 

+ 0.598 jr- X)]0 + [C T (C5)(0.639A + 0.598 *2. x) ] 6 

L L C c 

c a c 

c 


( 94 ) 


U 


U l 


V 


= -[C (C6) (0 . 639A + 0.598 = 2 -)]6 | . - [G 2 (A + 0.598 -p) (^A + 1) ]-^ 


c o 


Equation (94) is the second in the set of linearized equations 
describing the complete system. 

The inertial pitching moment about the trimmer hinge axis can be 
obtained from equation (41) of appendix A. Ignoring the lateral-directional 
variables, the nonlinear expression for the pitching motion of the trimmer 
is 

\ > I y „(« c + « p + q) + » c Xcg ( 1X ™ ( 'p + ’> + + ^ Sin 4 p 

aero 

+ [X CH ( % + c l )2 “ Z CH ( % + < ’ )]C0S 6 p " U - qW}sin(S p + 6 c ) 

? .. .. (95) 

+ m c X cg {[ - X CH ( 'p + 4 + Z- H ( 5p + q)]sin 6 p + [X* R (6 p + q) 

4- Z ^(6 + q) 2 ]cos 6^ - W 4* qU}cos(6^ 4- 6^) 

In linearized form, the fundamental trimmer pitching equation is, 

therefore 


aero 


M " [I „l« c + [In + m c P z iZ' cos « - sin « ) 

J J H o o 

c 

+ V*. (z ch sin { P + X CH cos 6 P )]6 P + [ V 
Ho o 

c 


(96) 


+ m„p„ (z; u cos - x' ru sin 6 p ) + i^p (Z^ R sin 6^ 

O n 


c k z tt ' CH 
h. 

C 


p CH 
r o 


+ X’ cos 6 )]q + [m p U ]q + [m p U ]U + [-m p U lot 

pll t-v * 1 /> r v r\ r* ^ n n -v o 


*H 


c z R o- 
c 


c 
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where 


p = X" cos (8 +8 ) 

*H Cg P o C o 

c 


p = -X" sin (8 +8 ) 


P c 
o o 


Aside from the circulatory lift of equation (94), a lift force arises 
from the apparent mass of air accelerated by the trimmer motion. This lift 
force, according to reference 5 , is 

pU s c 

L = ■ - a (98) 

m E c 

c 

Differentiating equation (83) , 

a = a +8 +8 + (Cl)8 + (C2)q 4- (C3)8 4- (C4)C (99) 

c F p c p c 

To compute the pitching moment acting on the trimmer, the 
apparent-mass lift is segregated into components arising from plunging accel- 
eration (acting at the one-half chord point) and from pure rotation (acting 
at the three-quarter chord point) . 

Accordingly, the lift acting at the one-half chord point is 

pU S C 

L = — (a - 8 - 8 ) (100) 

m^^2 E c p c 

and the lift acting at the three-quarter chord point is 

pus c 

L = (6 + 6 ) (101) 

m 3 /4 E P c 

The forces acting on the free surface are shown in figure 22. 


L-circ 


L m |/2 L m 3/4 



K q * Hinge axis (shown in a 

~~}T ** • position for clarity in 


position to: 
derivation) 


FIGURE 22. AERODYNAMIC FORCES ACTING ON FREE SURFACES 
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From the figure, the increment in aerodynamic pitching moment is: 


M 


A C A C 

= M„ 6 + L . (X + -f) + L X + L (X - 7-% 

S t 1 olrc c 4 ”1/2 c "3/4 c 4 

aero 


(102) 


or 


pU 2 _ pU 2 A c 

M y = t C m, -f S c^ 6 t + hr S c< X c + T* ]C L 

y c S c 

aero t 


+ t- 


pU s c 


A pU s c 

X ](a - 6 - 6 ) + [ °_ C - ■■ 

c c p c E 


(X. 


(103) 


-7r )1(6 p + V 

Equations (103) and (96) can be combined. Combining terms and 
converting to operational notation, the trimmer pitching moment equation is: 


[ (C22) a + (C30 - C31) X] 6 + [-(C30)A]ct + [-C29]C T 

C C .L 


+ [ (C27) A]a + [ (C23) A + (C30 - C31)A]<$ + [(C24)A 

F p 

+ (C25) A ] 0 + [ (C26) A]u = + [C28R 


(104) 


where 


C22 - I y „ 

C23 = I + m o (Z' cos 8 - X ' sin S ) 

} H *o *o 

c 

+ VV, (Z CH Sin *p +X CH 008 S P > 
n O O 

C 

C24 = C23 


C25 = m p U 

C Xg O 

c 

C26 = m p U 
c z H o 

c 

C27 = -C25 

C28 = C 1/2 pU 2 S C 
til. o c c 

6 t 


(105) 
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C29 = 1/2 pU 2 S (X + 
o c 



C30 

C31 


pU S C 
o c c 


pU s c 

o c c 


X 

c 



Equation (104) is the third equation in the set describing the 
complete system. 

The physical coupling between the free- trimmer and the wing is 

implicitly incorporated into equations (85), (94), and (104), but to depict 

the effect of the trimmer on the wing/boom assembly, the force acting at 

the trimmer hinge must be evaluated. 

From equation (75) , the force transmitted from the trimmer to the 

boom (F , ) can be calculated, since 
c->b ’ 


' c-*b 


= -F, = AF 

b-*c aero 


+ Am g - m a 
c c c 


(106) 


c 

If the aerodynamic force is computed in the wing hinge axis system, 
the trimmer lift force can be defined as acting solely in the -z^ direction and 
the drag force as acting solely in the -x^ direction. 

Since we are dealing with perturbations from the nominal values in 
the linearized equations. 


AF = -AL 1 - AD 1 

aero c z c x 

c H H 


(107) 


The weight force, expressed in the wing hinge axis system, is 


m g = [-m g sin 0]1 + [m g cos 0 ] 1 

c c *h c Z H 


So, in linearized form. 


(108) 


A (mg) = [-m g6]l 
C 


(109) 


From equations (109) , (107) , and (106) , the x-component of the force 
transmitted from the trimmer to the boom is 


F 

X c->b 


= -AD 


c 


m g0 - m (a ) 
c° c ccg X 


(110) 


The drag of the surface is 
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f 


dC D 2 2 

= [C n + ~ T C t n/2 PU S 

c D 2 L o c 

o dCL c 

c L 


Taking differentials, the linearized form is 

+ 


40 c - [ (C D 


dC D 2 2 ^C 2 

-2 C L )cU o S c'“ + t-l 
dC o dC T o c 

L c L c 


( 111 ) 


( 112 ) 


'where 


Incorporating equations (112) and (74) into equation (110), 
F 


c+b 


= [C7]u + [C9](q + 5 ) + [C10](q +6 + 6 ) + [Cll]u 

P P c 

+ [C12 ]C + [C13]e 

Li 

C 


C7 =-m U 
c o 

c9 = ^C (Z CH COS 6 p ‘ X CH Sin 6 p } 

o r o 

CIO =-m [-X" sin(<5 +6 ) + Z" cos(6 +6 )] 

c c § Po C o cg P o °o 

dC D 2 2 

cn.- [c +-jc ]pu o s 

o dC T o 

c L c 


dC 


C12 = - 


£• C t PU 2 S 
2 L o c 


(113) 


(114) 


dC T o 
L c 

C13 = -m c g 

Finally, in operational notation, the x-component of the force trans- 
mitted from the trimmer to the boom is 

[-1]F + [C12 ]C + [ (C9 + CIO) X 2 ] 6 + [C13 + (C9 + ClO)X 2 ]0 

x c+b L c P 


(115) 


+ [(C10)A ]<S + [Cll + (C7)A]u = 0 

c 


Equation (115) is the fourth of the set of equations describing 
the total system. 

Similarly, the z, component of the hinee force is 
h 

F = -AL - m (a ) 
z , c c ccg z 

c+b & 


(116) 


i 
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The increment of lift has contributions from the circulation lift 
(equation (92)), the apparent-mass lift (equation (98)), and airspeed 


changes : 


7 pU S C 9 

AL c - q [1/2 pD o S c ] + a c + [C L pU o S c ]u 

c o 

c 

Substituting equations (117) and (74) into equation (116) , 


(117) 


where 


F = [C14]d + [C15]q + [C17 ] (q + 6 ) + [C18](q + 6 + 6 ) (118) 


+ [C19]C + [C20]d + [C21]u 

JL C 

C 


C14 = -m U 
c o 

C15 = -C14 

C17 = m (Z* sin 6 + X' cos 6 ) 

c CH p CH p 

o *o 

C18 = m [X" cos (5 +6 )] 

c cg P o C o 


C19 = -1/2 pU S 
o c 


C20 = - 


pu 2 s c 

O C C 


(119) 


C21 = -C T pU S 
L o c 
o 

c 


or, in operational notation. 


[-1]F + [(C14)A]a + [ (C15 )A + (C17 + C18)X ]6 

Z c-b F 

+ [ (C17 + C18)X 2 ]<5 + [ (C18) X 2 ] 6 + [C19]C 

P CL 


( 120 ) 


+ [(C20)X]a c + [C21]u = 0 

Equation (120) is the fifth in the set of linear equations describing 
the total system. 


Wing/Boom Assembly Free Body 

The nonlinear expressions for the acceleration of the center of 
gravity of the wing/boom assembly are obtained by eliminating the lateral- 
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► 


irectional variables from equation (28) in appendix A. 

(a, ) = U + qW - (q + 6 ) 2 (X' cos 6 + Z' sin 6 ) 

WB x p eg p eg p 

eg 

+ (6 + q) (-X' sin <5 + Z' cos 6 ) 

P eg P eg p 

(a__ ) = W - qU - (q + <5 ) 2 (-X' sin <S + Z’ cos <5 ) 

WB z p eg p eg D 


(121a) 


eg 


- (6 + q) (X' cos <5 + V sin 6 ) 

P eg P eg p' 


(121b) 


In linearized form, these expressions become 

(a,— ) = U + [-X' sin 6 + Z' cos 6 ] (q + 6 ) 

WB x eg p eg P P 


eg 


(a^ TB ) = W + [-U ]q + [-X' cos 5 - Z' sin 6 ] (q + 6 ) 

WB__ z o eg p eg p M p 


( 122 ) 


eg - P o Cg V ‘ P ' 

Considering the wing/boom assembly as a free body, the motion is 
governed by the aerodynamic lift and drag forces, the force transmitted from 
the free trimmer to the boom through the trimmer hinge, the force transmitted 
from the fuselage to the wing through the wing hinge, and the force of 
gravity. Thus, 


F + F , + F. 

aer °WB 


c-b + F F->W + "V# "V^WB 


(123) 


eg 


as 


The aerodynamic force on the wing/boom assembly can be represented 

(124) 


F =“AD 1 - ALl 

aero WB x TJ WB 2 


"WB ““ X H 
By analogy to equation (112) , 

dC 


H 


dC T 


4 d wb" [(c d + ZT c t )P^S]u+ [-B-C l ptfW 
°w dC L „ °W dC L_, °W ^ 


( 125 ) 


J W - \ 

By analogy to equation (82) , the geometrical angle of attack of 
the wing, in linearized form, based upon the normal velocity at the one-half 
chord point, is 

“w = + 6 p + f COS 6 P ](q + V + ( dT~ )C L (126) 

r o o L c 

c 

The last term in equation (126) is the induced effective angle of 
attack increment caused by the vortex system on the free trimmer system. 
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Equation (126) will be written as 




= 


+ 6 + [G33 ] (q + 6 ) + [C34]C 


(127 


where 


C33 = - — cos 6 
U p 

o r o 


da 


C34 = 


W 


dC 


(128) 


So, in operational notation, 

[l]a w + [l]a F + [1 + (C33) X ] 6 + [(C33)X]0 + [C34]C L = 0 (129), 

P c 

Equation (129) is the sixth in the set of linearized equations i 

representing the system. 

As with the trimmer surface, the circulatory lift coefficient of ! 

the wing has contributions from the motion of the airplane and from the 1 

vertical gust velocity. 

V 

C L = G 1 (t)[o w + (C36) (q + & )] + G 2 (t) ^ (130) 


where 


W 


w 


“W 


C36 = 


ttU 


(131) 


In operational notation, 


U 


0.639A + 0.598 ~ 

G 1 = C L t 

a W X + 0.598 


„ _ ^ r -, 0.4881 0.272X 

z " C L 11 

a W X + 0.455 X + 1.04 


0.193X 


U 


-] 


x + 4 - 71 r 


(132) 


(133) 


U U 

f x + 1 -°^s 

Equations (132) and (133) are the operational equivalents of the 
indicial response function given in reference 5 for a wing of aspect ratio 6. 

For computer mechanization of the homogeneous equations, it is 
desirable to eliminate the denominators of the coefficients. This is done by 
multiplying equation (130) by the denominator of equation (132) . 

The final expression for the wing lift coefficient is 
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(134) 


U U 

-[X + 0.598 + [C L (0.639A + 0.598 )]o y 


<*W 


u 


A(0.639A 


+ [(C36)C t A(0.639A + 0.598 =^)]<$ +[(C36)C T 

V C P \ 

“w a w 

U U V 

+ 0.598 — )]0 = t-(X + 0.598 ~)G 2 

Wo 

Equation (134) is the seventh in the set of linear equations 
describing the system. 

The total wing lift force is made up of the circulatory lift 
defined by the lift coefficient of equation (134), and the lift due to the 
acceleration of the apparent mass of air surrounding the wing. The apparent 
mass lift is 




pU sc 

o 


m 


E “W 


pU SC 

-f- <*w - y + 

Acts at 1/2 
chord 


pU SC . 

— § — 6 
E p 

Acts at 3/4 
chord 


(135) 


In linear form, the increment of lift on the wing is 


pU SC 


AL W = [ 1/2 PU 0 S]C + + [C L P%S]u 

^ °W 


(136) 


Equations (136), (125), (115), (120), and (122) can be substituted into 
equation (123) to obtain explicit equations for the components of the force 
transmitted from the wing to the fuselage through the wing hinge. In opera- 
tional notation, these are: 




[-1]F + [C38 + (C40) A]u + [C39]C T + [1]F 

Vf 

+ [C43 + (C42)X 2 ]0 + [(C42)X 2 ]<S p = 0 


c->b 


(137) 


[~1]F + [C44]C + [ (C45) X]cl + [C46]u 

Vf h ^ 


+ [ 1 ] F + [(C47)A]<x_ + [ (-C47) X + (C49)X ]0 

Vb F 


(138) 


+ [ (C49) A ] 6 = 0 

P 
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where 


D 2 2 

C38 =-(C + C ) P U s 

U L o 

% d \, % 


D 2 

C39 7T- C T pU S 

Ar^ 0 

\ « 


040 ■ -”WB°0 

C42 ’ -“wB^^g Sl " 6 p + Z cg COS 5 p > 


C43 = “i 


C44 = - 1/2 pU s 
o 


C45 = 


pU SC 
o 


C46 = -C pU S 
L o 

°W 


C47 = C40 


C49 = m (X T cos 6 4- Z f sin 6 ) 

WB eg p eg p , 


Wi5 eg P o eg p Q ', 

Equations (137) and (138) are the eighth and ninth equations in the 
linear set describing the system. 

The pitching moment about the wing hinge caused by the force 
applied to the boom by the free trimmer is the cross-product of the vector 
from the wing hinge to the trimmer hinge and the applied force. The magni- 
tude is 

“f ■ ( - X iH Sln S p + Z i« COS 5 p > F x v 
c r o *o c->b 

- < X i H 5 p + Z i H Sl ” S p > F Z „ <UC 

o *o c+b 

For small displacements, the increment in pitching moment about 
the hinge axis due to gravity is 


ight ■ + y 
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| By analogy to equation (103) , but allowing for the fact that the 

wing boom aerodynamic center may not be at the quarter-chord point, the 
aerodynamic moment on the wing is 


M 


— X 

. C ac Mm , r o 


pU SC 


aero 


- [1/2 + [—£— X1 (5 W - y 


PU SC _ x . 

+ t-f— « - 1 6 p 


(142) 


The inertial pitching moment on the wing boom assembly is obtained 
from equation (47) in appendix A by eliminating the lateral-directional terms: 


M 


where 


^inertial 


= I ,(6 + q) - iil_[(W - qU)(X' cos 6 + Z' sin 6 ) 

y P wB eg p eg V 


eg 


eg 


- (U + qW) (-X* sin 6 + Z' cos 6 )] 

eg p eg p' 


In linearized form, this becomes: 


M 


^inertial 


- I y .« p + “o'’ + Vz V* 


y p 

- [m WB P * °o 14 F 
eg 


Cg 


eg 


p = X cos 6 + Z sxn 6 

x eg p eg p 

eg *o 65 *o 


p = -X f sin 6 + Z f cos <5 

Z cg Cg P o cg P o 


(143) 


(144) 


(145) 


Combining equations (140) , (141) , (142) , and (143) , the pitching 
moment equation for the wing boom assembly is, in operational form: 

[C51 + (C56 - C59) A + (C50)A 2 ]<5 + [C51 + (C52)A + (C56)A 2 ]8 

p 


+ [ (C53) A]u + [C54]C + f- (C52) A]ot_ + [(C55)A]a 

L r W 

W 


(146) 


+ [C5 7 ] F + [C58]F = 0 

X c-»b Z c->b 

Equation (146) is the tenth in the series of linear equations where 


C50 = I , 

y 


C51 = 


m WB gZ c 


eg 


C52 ■ W* D o 

Cg 


(147a) 


103 


-1/2 pU S(X + -=■ 
o l 

pU SC „ 

o 



o 

E 

X > 

1 

f) 

-1 

sin 

6 

- Z* 

L CH 

P o 

CH 

- ? 

cos 

6 

+ Z' 

L CH 

p~ 

Cl 
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Fuselage Assembly Free Body 

The inertial pitching moment about the wing hinge axis is given 
by equation (54) in appendix A. For longitudinal motion, the expression may 
be written 

M y = I q + m F [(U+ qW)Z p - (W - qU)^] (14! 

^inertial ^ 

In linear form, 

M = I q + [y F U 0 ^ + [-Wo^F + [ Vp U o ]q (14 

Inertial F 

The applied moments on the fuselage assembly arise from aerodynamic 
moments and the weight moment caused by the offset center of gravity. 

The weight force vector is 

F p = (-mpg sin 0)1 + (mpg cos 0) (15 

rl H 

The pitching moment caused by the weight force is the cross-product 
of the position vector from the wing hinge to the fuselage center of gravity 
and the weight force vector: 

= _m pg z F sin 0 - C ° S 6 ^ 15 

^weight 

Or, in linearized form, with 0 =0, 



M 


[ ^g^p] 


0 


(152) 


weight 

The aerodynamic moments on the fuselage will be assumed to have con- 
tributions from ot„, q, 6 , V , and V . 

* P g g 


where 


For the angle of attack contribution, 

AM = C 1/2 pU 2 SC(a_ + a ) 

aero m o F gust 

a « F 


V JL V 
= _£ _ _£ 


gust U 


u u 

o o 


(153) 


(154) 


Since most of the moment arises from the horizontal tail on the 
i fuselage, it is reasonable to use the gust velocity adjusted for the penetra- 
tion time from the wing to the tail. Therefore, the pitching moment contribu- 
tion from fuselage angle of attack is 

2 - 


2 - V Z t V 

AM = 1/2 pU SCC (a_ + -£■ - -£• ) 

aero o m F U U U 

a o o o 

F 


(155) 


The aerodynamic pitch damping term includes a contribution from 


the gust gradient 


AM = 1/4 pU SCC (q - 77 ^) 

aero o m n U 

q q ° 


(156) 


The downwash angle at the horizontal tail can be expressed in terms 
of wing circulatory lift coefficient as follows: 

de 1 


e = 


da C T 


(157) 


w 


Thus, the effect of wing lift on the fuselage pitching moment will 


( be 


de 

AM = [C 1/2 pU 2 SC 7 r^-]C T 

aero c m £ o C L L y 

L a W 


(158) 


Combining the applied moments and setting those equal to the iner- 
tial pitching moment, the fuselage pitching equation becomes, in operational 
notation. 
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[C 62 + (C65 - C61) X - (C59)X ]0 + [(C60)X]u + [C63 + (C61)X]a. 


+ [C66 ]C + { [ C63] +[(C64)(C63) - C65]A}- s -= 0 
L W o 


C59 = I 


C60 = -m F Z p U o 


C61 = m. 


pVo 


C62 = -m F gZ F 

C63 = 1/2 pU 2 SCC 
o m 


C64 = - 


C65 = 1/4 pU S(C) 2 C 
o m 


C 66 = (C 1/2 pU 2 SC) ( 7 ^-) 
m o C_ 

e L 

a w 


Equation (159) is the eleventh in the set of linear equations des- 
cribing the system. 

The components of the acceleration of the center of gravity of the 
fuselage can be expressed in the wing hinge axis system as: 

(a p ) = U + qW + Z p q - q 2 X p 

Cg ( 1 ( 

(a F ) z = W - qU - Xpq - q 2 Z p 
Cg 

In linearized form, with q = w =0. 

n o o 

(a F ) x - U + Z F < 


(a F ) 2 ' W - U Q q - J^q 
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The dynamic equation for motion in the direction is. 


V a F = F x 

^ F cg 


+ AF 


aero 


mpge 


And, 


AF x = -AD p = -[C D pU Q S]u 
aero F 


(163) 


(164) 


Substituting equations (164) and (162) into (163) , the airspeed 
equation is obtained. In operational form, this is: 

[C69 - (C67) X]u + [C70 + (C68)X 2 ]6 + [1]F = 0 

X W>F 


(165) 


where 


C67 = m^ 
C68 = -mpZp 


(166) 


C69 = ■ -C pU S 
D F ° 

C70 = -n^g 

Equation (165) is the twelfth in the set of equations describing 
the total system. 

The equation of motion in the direction is 


”F <a F >2 " V„ + 


AF 

eg Z aero 


167) 


The aerodynamic force increment in the direction is assumed to 


be dependent upon a^, u, and V^. 


AF„ 


v 

= - [C T 1/2 pU 2 S]a F + [C T pU 2 S]u + [C T 1/2 pU^S]^ (168) 


aero 


l f ■ ° 

O 


o J U 


"F "o F 

Substituting equation (168) and (162) into (167) , and writing in 
operational form, the thirteenth and last of the set of linear equations is 
obtained: 


[C72 - (C67) A]a + [(C67)A + (C71)A ]0 + [1]F 

r A 


W+F 


V 

+ [C73]u = -[C72] ^ 

o 


(169) 


where 
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C72 = -C T 1/2 pU S 
L o 

“F 

C73 = -C T pU 2 S 
L f o 

O 


Matrix Form 

As derived, the state vector of the system has 13 components 


: F 


x 


c->b 


J c+b 


W 


X = C 


W 


V»-F 


"W-KF 


0 

u 

a. 


I 


_ F j 

In matrix form, the equations of motion are: 

V 


[A]X = [B] 6^ + [C] ^ 



In the last equation, [A] is a 13 x 13 matrix of the coefficients 
of the homogeneous equations, and [B] and [C] are column matrices through 
which the control tab displacement and gust velocities are introduced to 
perturb the system. 

Conversion to First-Order Form 


The numerical integration of the equations of motion requires that 
an explicit solution be developed for the highest derivative required for 
| each variable. 

Although two of the variables of equation (172) , and a^, 

I appear in nonderivative form, their first derivatives are used elsewhere 

| in the set of equations because of the unsteady aerodynamic functions. 

To modify equation (172) for direct numerical integration, the 

a and a T _ equations were replaced by their time derivatives, and all lower 
c W 

order derivative terms were moved to the right side of the equality sign. 

In the computer subroutine for evaluating the derivatives at each 
instant of time, it was necessary to perform a simultaneous solution for 
the 13 highest-order derivatives of the modified equations. This solution, 
combined with another 12 direct integrations for the lower order variables, 
yielded the current values of all quantities needed in the integration 
process. 

Aerodynamic Coefficients 

| The aerodynamic coefficients of the representative aircraft were 

j estimated from two primary sources. A vortex-lattice program, modified at 
I Battelle Columbus Laboratories by using the method described in reference 7, 

i 

| was used to predict the static wing and free-trimmer characteristics; the 
approach outlined in reference 8 was used to estimate the primary contribu- 
tions of the fuselage and fuselage-mounted stabilizing surface. 

Geometrical data for the analyses were obtained from the three- 
view drawings of the candidate aircraft contained in the main body of the 
report . 


109 



Static Wing/Trimmer Coefficients 


Lift-curve slope . The lift-curve slopes for the isolated wing and 
trimmer surfaces were obtained from the vortex lattice program by using 12 
spanwise vortex elements per semispan. For aspect ratio 6, the resulting 
value was 4.276 per radian. 

Profile drag coefficients . The section profile drag coefficient 
for both surfaces was broadly estimated to be 0.01. 


Mutual interference coefficients . Because of their close proximity, 
the wing and trimmer have significant mutual interference effects. To account 
for this phenomenon, the effective induced angle of attack on each surface 
was computed as a function of the lift coefficient on the other surface. 

Using the appropriate relative geometry of the surfaces, the 
influence coefficients were computed from: 


da 

c 




1 



da 


W 


dC 


L 


c 





(173) 


In equations (173), the prime superscript denotes the value obtained 
without the presence of the other lifting surface. 

Representative values, so obtained, are given in table XV for an 
area ratio of 1/6. 


TABLE XV. MUTUAL INTERFERENCE COEFFICIENTS 





da 

dC 


W 

L 


c 


2.0 -0.50 0.0226 -0.0140 

0.75 -0.364 0.0778 -0.0122 

-0.75* 0 0.0898 0.00585 

-1.50* 0 0.02772 0.004473 


* Tip-mounted aft trimmers. 
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I Induced drag coefficients . The induced drag coefficients for both 

[the wing and trimmer surfaces were estimated to be 0.0624. 

i 

Fuselage Coefficients 

The aerodynamic coefficients of the fuselage and fuselage-mounted 
horizontal stabilizer were estimated by assuming an equivalent circular body 
with area distribution as used in the aircraft shown in figure 2 of the main 
body of the report . 

Following the method of reference 8, the lift and moment contribu- 
tions of the fuselage and horizontal tail were determined. The reference 
point for the moment coefficients was the wing hinge axis. 

, Fuselage-tail lift-curve slope . The lift coefficient of the 

fuselage- tail assembly, as a function of fuselage angle of attack is estimated 
to be: 

C = 0.00751(a - 4°) + 0.000033(a^ -4°) + (0 . 0104) (0 . 55aJ (174) 

L F f f f 

From which C =0.758 per radian. 

L 


Fuselage static angle of attack stability . Similarly, the slope of 
the fuselage moment coefficient with respect to angle of attack is 


! 


I 

I 

f 

i 



= -0.006 - 


0.00015a per degree 

r 


From which C = -0.344 per radian, 
m 


(175) 

(176) 


Pitch damping coefficient . The pitch damping coefficient was esti- 
mated on the assumption that all damping arises from the forces on the 

horizontal tail. On this basis, C = -4.76 per radian. 

m 


111 



APPENDIX C. METHOD OF COMPUTING TURBULENCE RESPONSES 


x 


Symbols 

scale length of turbulence, meters (feet) 
airspeed , meters/second (feet/ second) 
rms value of variable x 
power spectral density function 
spatial frequency, radians/meter 


Responses to Continuous Turbulence 


1 


Equation (1) in this report describes the deterministic response of | 
the longitudinal system to the vertical gust velocity. 

For random turbulence responses, the frequency response function is 
used to compute the spectrum of the response for each variable of interest. The 
output spectrum for a variable, x, is given by the following equation (ref 6): 


$ x (fi) = 


x 


g 


*(n) 


(177) 


x 


where — is the modulus of a frequency response function which defines the re- 
v g 

sponse of the variable to the gust velocity. 

The root-mean-square response of the variable is the quantity of 
interest, and it is then computed from: 

1 1/ 2 


a x = 


/ $ v (fl) dQ 

O x 


(173) 


In the numerical integration of equation (178) , the actual limits of 
integration were from spatial frequencies corresponding to temporal frequencies 
ranging from 0.3 to 40 radians per second, using the relationship 


a = 


u 


(179) 


Only vertical gust components were considered, and the Dryden power 
spectral density function was used with a scale length of 533.4 meters 
(1750 feet) and an rms gust intensity of 0.305 meters per second (1 foot per 
second) . The Dryden spectrum is given by 


*(fl) = 0.0283 0 2 - 1 + 3 ^ 2L - 2 

8 * (l+fi2 L 2)2 


(180) 
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The frequency response function was computed directly at each fre- 

I 

quency by setting the gust velocity equal to j oo in the equations of motion. 

In this way, the steady-state sinusoidal response of each variable of interest 
was computed for each value of the sinusoidal gust velocity frequency, u). The 
absolute magnitude of the response function was then computed, as a function 
of frequency, to obtain the required frequency response function for use in 
equation (177). 
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